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ABSTRACT 

This work investigates the effects of digital processing 
in radiometers. It deals mainly with two digital versions 
of a total power radiometer. The first consists of RP, 

Mixer and IP sections followed by an analog to digital 
converter. All further processing is done in a digital 
computer. The second version consists of RP, Mixer and IP 
sections followed by a square lav; detector, RC filter and 
analog to digital converter. Prom this point on the pro- 
cessing is done by a digital computer. A figure of merit 
is defined based on the performance of an analog total power 
radiometer. Exact results are obtained for the figure of 
merit of the first digital version. Por the second, an 
approximate solution is obtained. The effects of saturation 
and finite step size of the quantizer were taken into con- 
sideration for the above results. The performance of digital 
balanced-Dlcke and noise-adding radiometers is investigated 
using, the above results. The effects of digital filtering 

V 

on the performance of a radiometer is considered. ' 
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I. INTRODUCTION 



A radiometer is basically a highly sensitive and stable 
noise receiver. The principle of operation lies In the fact 
that any object at temperatures above absolute zero radiates 
energy In the form of electromagnetic waves. Until some 
years ago, radiometers were used mainly as radio telescope 
receivers and constituted the main tool of Radio Astronomy 
[Ref. 1]. Lately, their use has expanded to many other 
areas. The measurement of atmospheric temperature [Refs. 2 
and 3]> detection of air turbulence In the Troposphere 
[Ref. 4], airborne mapping [Ref. 5], the measurement of 
absolute radiation from a projectile flow field [Ref. 6], 
and passive detection [Ref. 7] can be cited as examples. 

The type of processing Involved In a radiometer Is 
especially .fitted to digital methods. A number of Insti- 
tutions have Integrated the radiometer with the digital 
computer In their applications . The purpose of this work 
Is to Investigate the effects of this Integration, and to 
provide guidelines to the design of digital radiometers. 

The main body of the Investigation will deal with two 
digital versions of a total power radiometer (TPR) shown In 
Figure 1. The first one (Figure 2) consists of RF, Mixer, 
and IF sections, followed by an Analog to Digital Converter 
(ADC) section. All further processing is done In a digital 
computer. The second one (Figure 3) consists of RF, Mixer, 
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and IF sections followed by a square law detector, a low 
pass filter, and an ADC. From that point on, the processing 
is done by a digital computer. 

The main difference of the two cases is the location of 

I 

the ADC. The second version (Figure 3 ) is presently being 
used and its main advantages over the first version (Figure 
2) is that it requires a much lower sampling rate. 

The results obtained for a total power digital radiometer 
can be used to determine performance of other types of 
radiometers. Performance of digital Dicke and noise-adding 
radiometers will be considered. 







FIGURE 1 
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II. PERFORMANCE FIGURE OF A RADIOMETER. 
DEGRADATION FACTOR 



By definition, the sensitivity, or minimum detectable 
temperature, of a radiometer is taken as that change in 
operating temperature T that causes a change in the expected 
value of the output equal to one standard deviation of the 
output [Ref. 1]. Refer to Figure 4. It shows the expected 
value and the standard deviation of the output voltage as 
functions of temperature for a typical radiometer. 




FIGURE 4 
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From Figure 4, It follows that (to first order) 



AE[I] = 



dECi] 

dT 




AT 



( 1 ) 



where E[*] = Expected value of [•], and T^p Is operating 
temperature. But by definition 



Therefore 



AE[I] = OjL 

' op 



= dE[I] 
dT 




( 2 ) 



AT = 



dETTF T 
dT op 



(3) 



A performance figure can be defined as. 




0 , 





(4) 



It can be shown [Ref. 1] that for an analog total power 
radiometer (Figure 1), using the fact that the Inverse of 
the Integration time Is much smaller than the bandwidth of 
the RF, Mixer, IF filter, that 



E[I] = TkT^p / lH^(f)l^ df (5) 

OO I 

( 5 ^ = ^ / |H^(f)r df (6) 
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where 



T = Integration time 



H^(f) 



Band-pass filter (RP,Mlx,IF) frequency 
response 



k = Boltzmann's constant 



Substituting into (4) 



AT 

T 

op 






(7) 



where 






( / af)^ 

/ |H (f)|!* af 
0 



( 8 ) 



The performance figure of the analog TPR, (7), will be 
used as a basis for evaluating the performance of a digital 
TPR. Hence, the Degradation Factor can be defined as 



p _ Performance figure of Digital Radiometer 
Performance figure of Analog Radiometer 

Intuitively, F should be equal or greater than one since 
as the sampling rate of the ADC goes to infinity and the 
step size of the quantizer goes to zero, both digital 
radiometers are reduced to an analog TPR. 

Some comments are appropriate at this point. First, 

(5) and (6) show that for analog radiometers E[I] and 
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are linear functions of temperature and go through the 
origin. Hence equation (^) reduces to 



AT _ "I 

t—'eTTI 

op 

The addition of a quantizer, due to saturation, makes 
E[I] and nonlinear with respect to temperature, and the 
performance figure is thus temperature dependent. There- 
fore (4) must be used instead of the better-known (7). 
Secondly, the parameter T^^ in equation (^) is not a 
convenient choice, since for a given T^^, different gain 
settings of the RF, Mixer, IF filter will change the 
performance of the radiometer. Referring to Figures 2 and 
3 it is easily shown that 



■ • “ "x 

c 


(Figure 


2) 


(10) 


T = S 0^ 


(Figure 


3) 


(11) 


where a and 3 are proportionality 


constants . 


Using the 



2 

parameters cr^ or cr^ as appropriate, instead of T, norma- 
lizes the problem of gain setting out of the performance 
figure . 
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Then, for the radiometer of Figure 2 




dELU 

da^ ^ 



^ at T 
c 



op 



( 12 ) 



and for the radiometer of Figure 3 



AT 



op 



dELl] 



at 



T . 

op 



(13) 
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III. DEGRADATION FACTOR CALCULATION WITH NO QUANTIZER 



An ADC performs two basic operations on a signal. One 
is sampling, and the other is amplitude quantization. In 
this chapter the effects of sampling alone are considered. 
Hence it is assumed that each sample can be represented 
without any error. 

\ 

A. IP SAMPLING 

Refer to Figure 2. Assume that at the antenna white 
gaussian noise is present with double-sided power spectral 
density No/2 watts/hertz. Then X(t) is a Narrow Band 
Gaussian Process with known power spectral density, and 
can be represented as 

X(t) = X (t)cos2TTf t - X (t)sin2iTf t (l4) 

.c o s o 

where X (t) and X (t) are baseband Gaussian processes and 
s c 

f^ is the IF center frequency of the radiometer. It can 

be shown [Ref. 8] that if X(t), X (t) and X (t) are to be 

c s 

wide sense stationary, then 

Rx (y) = Rx (15) 

■^c "^s 

Ry(p) = Ry (y)cos2iTf y + Ry Y (y)sin2iTf u (l6) 

C SC 
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where R (u) is the autocorrelation function of the process 

X(t) and Ry ^ (u) is the crosscorrelation function of the 
■^s c 

'processes X (t) and X (t). 

c s 

Furthermore if H (f) is symmetrical about f (this is 

o o 

the only case considered in this work) , then 



Rx X 
c s 



(17) 



and R„ (y) is given by 
^c 

oo 

R^ (y) = 2 / S^(f)cos[27r(f-f )y]df (l8) 

■^c 0 ^ ° 

and S„ (f), the power spectral density of X (t), is the 

A C 

C 

Fourier Transform of equation (l8). 

Physically, for the case of Interest (l.e. X(t), X^(t) 

and X (t) wide sense stationary, and H (f) symmetrical about 
s o 

f ), S„ (f) is Just twice the low pass version of S„(f). 

O A A 

C 

Referring to Figure 2, the low pass filters filter out the 
double-frequency terms coming out of the multipliers. Then 

\ 

V 

I = + I 2 (19) 

I, = ^ r X^^(kA) (20) 

1 n k=i c 

since X (kA) = q (kA) and 

c c 
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( 21 ) 



n 2 

^2 n ^ x/(kA) 

k=i s 

since X (kA) = q^(kA). 
s s 



Using the fact that for Gaussian random variables 



E[X^X2X^X^] = E[X^X2]E0C^X^]+E[X^X^]E[X2X^]+E[X^X^]E[X2X^] 

( 22 ) 

and expanding the double summation, it follows that 

4 



2o 



X^ h n-1 . 5 

+ - Z (l--)R^ (kA) 
n n , T n X ^ 

k=l c 



( 23 ) 



and 



E[Ii] = <1^ 



(210. 



For the random variable I 2 the same results are obtained 

since (p) = R^ (u) 
c s 



Then 



E[I']= 20^ 



( 25 ) 



and 



Oj^ = + 2E[I^l2] - 2E[I^]E[l2] 



( 26 ) 
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Expanding E[I^l 2 ] and using the fact that R„ „ (p) = 0, 

■^c s 

the last two terms of ( 26 ) cancel out. 

Using equations (25) and (26) In (12) It follows that 



(_^)2 = 1 + 
op 



2 ^ ^ k 2 

Z (1 - |)R/(kA) 

na^^ k=l ^ c 



(27) 



or 



(_^)2 = A ^ 
op 






T 



TOv k=l 



( 28 ) 



where 



T = nA 



(29) 



The degradation factor squared Is 



„ 2B„A n-1 . . „ 

= BgA + — ^ Z (1 - i^)Rj^‘^(kA) 

Ov k=l c 



(30) 



Ry (m) Is a function of H (f), therefore (30) will be used 
with four different frequency characteristics. 

1. Ideal Bandpass Filter 
Refer to Figure 5 



N B slniTBp 



( 31 ) 



20 
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and 



B (32) 

» 

Substituting (31) and (32) into (30) 

= B„A + 2B„A Z^(l— ) ( 33 ) 

^ ^ k=l (TTkAB)^ 

But 

Bg = B = W (34) 

where Bg Is the equivalent bandwidth from ( 8 ) and W is the 
half- power bandwidth. 

Then (33) can be rewritten as 

p 2 = VL + 5 ln!lkW . 2 (WA)^, kslnfjkW ( 35 , 

k=l (irkWA)^ k=l (irkV/A)"^ 

Since t>>A, it can be shown that the third term of (35) 
is much smaller than the second and thus can be Ignored. 
Letting A approach zero (sampling infinitely fast), (35) 
becomes 

= 2 W / dt ( 36 ) 

0 ( irWt ) 

and then since t Is large, the Integral's upper limit can 
be replaced by infinity, and the integral reduces to one, 
giving the correct result. 
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Actually, by Including the third term as A goes to zero, 
P will be slightly less than one. The reason for this is 
that in calculating the performance figure of an analog TPR 
in (7) the fact that ^ << B was used to get (5) and (6). 

Had this fact not been used, as A approached zero, P would 
approach one exactly. 

P will be plotted as a function of WA. Hence, for this 
case the Nyquist sampling rate (twice the highest frequency) 
corresponds to WA = 1 . 

The analysis done in this chapter is just a special case 
of the analysis that will be done in Chapter V, hence all 
results will be plotted in the figures of that chapter. The 
results obtained from (35) for an ideal bandpass filter are 
plotted In Plgures l8 through 23. It is of Interest to note 
that for W less than one (sampling faster than the Nyquist 
rate), P Is exactly equal to one, an intuitively satisfying 
result . 

2 . Single Tuned Pilter 

Por a Single Tuned Pilter 



IH^(f) 



1 + 4 (. 



f+f 



0\2 



W 



1 + 4(- 



f-f 



W 



0)2 



(37) 



where W is the half power bandwidth and W is much less than 

f . It follows that 
o 



Sy^if) = ^ lH^(f)l 



2 



(38) 
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and 



and 



N 

Sy (f) = ° 



1 + 



(39) 



From (39), using the Inverse Fourier Transform 



N W I 

(m) = 4- 



(40) 



and 



AT ^ A ^ 2A kAv -2-rrkWA 

% ■ ^ ^ k=l ■ 



(41) 



Equation (4l) can be put in closed form by letting 



A = e 



-2ttWA 



(42) 



Using the fact that [Ref. 9] 






k=l 



1 - A 



(43) 



.k _ A -nA^ + (n-l)A^'*'^ 



S kA*^ = 
k=l 



(1-A)' 



(44) 



and after some tedius but straightforward manipulations, it 
follows that 



(ffi^)^ = ^ (coth ttWA) - (cosech^ ■rrWA)(l - (45) 

op ^ ^ 

and 

2 

F = ABg(coth ttWA) - (cosech^ ttWA)(1 - (46) 

2tBe 
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It can be shovm that the equivalent bandwidth of a 
single tuned fulter is [Ref. 1] 

Bg = ttW (47) 

I 

where W is the half-power bandwidth. Furthermore, since 
t>>A, the second term is negligible with respect to the 
first, hence 

F^ = iTWA(coth ttWA) (48) 



To compare with the analog case, let A go to zero. 
Then F would become one, which is the correct result. 
Figures 24 through 29 show the results obtained for this 
filter. 

3. Gaussian Filter 

For a Gaussian Filter 



[H^(f)]^ 



-v( 

e 



f-f 

~w 



^) 



2 



+ 



v( 

e 



£■+^ 0,2 
w ^ 



(49) 



where v = 4 In 2 = 2.773» W is the half-power bandwidth, and 
W Is much less than f^. 

It follows that 



Sy.(f) 

■^c 



N e 
o 






(50) 



and 






e 



( 51 ) 
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Replacing equation (51) into (30) and using the fact 



that for this filter 



Bg = = 1.505 W ( 52 ) 



it follows that 



2 kA 

= I. 505 WA + 3 .OIWA Z (1-^) e ^ (53) 

k=l ^ 



Using the same argument as before (A<<t), (53) becomes 

„ n-1 -|(TTkWA)^ 

= I. 505 WA + 3 .OIW Z e (54) 

k=l 

To compare to the analog case let A approach zero. Then 
(54) becomes 

P T -|(TTkWA)2 

F = 3 .OIW / e (55) 

0 

and since t is large, the integral’s upper limit can be 
replaced by infinity, and the integral reduces to one, as 
it should. The numerical results are plotted in Figures 30 
through 35 . 

4 . Butterworth Filter ' 

The last filter to be considered is a Butterworth 
Filter of order two. For this case 



= b-T '56) 

where W is the half-power bandwidth, and W is much less than f^. 
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Then 



S„ (f) 



N 

1 + 16 ( 1 )'' 



(57) 



an(i 



V“ -I'-l 7 ^, ,UUw 

^ ( 2 ) 



Ry (y) 






(58) 



For this filter. It can be shown that 
(o\^ 

ttW = 1.48W 

® 3 



(59) 



Substituting ( 58 ) and (59) Into (30) It follows, 
after some effort that 



n— 1 ^ 

= l.ilSWA + 2 . 96 WA E ""^^[1 + sln(k7r(2)^WA] 

k=l 

_ 2^96^AW)f_ + sln(k7T(2)^WA] ( 60 ) 

k=l 



Again It can be argued that the third term Is negligible, 
therefore 



n-1 



■ hr 



= l.i|8WA + 2 . 96 WA E *^^[1 + sln(k7T(2)^WA)] ( 6 I) 

k=l 



Comparing to the analog case by letting A approach 
zero, ( 61 ) becomes 

h. 



X ^ 

P^ = 2 . 96 WA / *^^[1 + sln(7T(2)*^Wt)]dt 

0 



( 62 ) 
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Since T Is large, the integral’s upper limit can be replaced 
by infinity, and it reduces as it should to one. 

The numerical results for this filter are plotted 
in Figures 36 through 4l. 



B. RC FILTERING AND THEN SAMPLING 

Refer to Figure 3 . This type of digital TPR was studied 
by Ohlson [Ref. 10] for the case where no quantizer is 
present. Only the results and its implications are presented 
here. 

It can be shown that 



N 

E[Z] = g ^ 



(63) 



and 

(SI) 

% 

where g is a system gain parameter, N^ is as defined on 
page 17 , and tj^^ is the time constant of the RC filter. 
Also 



E[I] = E[Z] 



(65) 



and 





n-1 

Z (1 - 
k=l 



^)[R2(kA)-E^[Z]] 



( 66 ) 



where ^^^(u) is theautocorrelatlon function of the process Z(t). 
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Substituting (65) and (66) into (13), we have 



op 



n-1 



kA 






r 1 _L 2A ^ kA N top -1 

[ - + — E (1 - — ) e RC ] 
" ^ k=l ^ 



. ( 67 ) 



Defining 






( 68 ) 



and using equations (^3) and (44) it can be shown that 



ir.2 A *.1, A 

F = Trr coth 



( 1 -e"'^/'*^RC)(_^ csch oi^)^ 



2t 



2tRG 



RC 



2t 



2tRC 



RC 






(69) 



Equation ( 69 ) is of the same form as equation (47). 

Note that in that case it was argued that the second term 

0 

was negligible since tW was large (oh the order of 10°). 

Here the same assumption can not be made since r-^ << W 

tRC 

(i.e. the RC filter bandwidth is much smaller than the 
bandwidth of the RF, Mixer, IF filter). The first term of 
equation (69) reduces to one as A approaches zero. The 
second term however, contributes negatively. Therefore F 
can be less .than one. The reason for this lies in the addi- 
tion of an RC filter to the radiometer,’ which has the effect 
of lengthening the integration time t. As A approaches zero, 
a comparison is being made of two TPR, one with an RC filter 
after the square law detector and the other without one. 

The one with the RC filter will have greater sensitivity 
(smaller AT) due to a longer effective integration time. 
Therefore F can be less than one. 
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However, if the effects of sampling only are of Interest 
(not the effect of the addition of the RC filter), the 
second term must be disregarded by assuming that f/tj^Q = 

This implies that all the filtering is done by the integrator 
and none by the RC filter. 

As explained in part A of this chapter, the results are 
presented in Chapter V. Refer to Figures 46 through 51. 

In summary, the effects of an ideal ADC (i.e. no 
saturation and step size equal to zero) was studied in this 
chapter. Its main importance is the fact that it provides 
a lower bound for the Degradation Factor and, that under 
certain conditions, a correction can be added to these 
results to account for the finite step size of the quantizer. 



V 
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IV. TRANSFER CHARACTERISTICS OF A QUANTIZER 
WITH A ZERO MEAN GAUSSIAN INPUT PROCESS 



In this chapter the effects of a quantizer on a zero 
mean Gaussian Input process are analyzed. Two cases are 
Investigated. Refer to Figure 6. It Is desired to find 
the autocorrelation function at point 2 given the auto- 
correlation function of X(t) and the knowledge that It Is 
Gaussian and has zero mean. Several approaches have been 
explored for the case of a quantizer alone. The one sug- 
gested by Kellog [Ref. 11] was used since It can be easily 
extended for the case of a quantizer followed by a square 
law detector, and also the equations that It leads to are 
In a form suitable for computer programming. Figure 7 
shows two quantizers, one with even number of levels and 
the other with odd number of levels. 

Let N be the number of levels of the quantizer. Then 



X 



1 



2L(1 - N/2) 
N-1 



(70) 



and 



^ 2L[1 - (N-H)/2] 
N-1 



(71) 
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FIGURE 6 



A. QUANTIZER ALONE 

The autocorrelation function of q(t) Is 



Rq(y) = E{Q[X(t)]Q[X(t+y]} 



Let 



X(t) = 



and 



X(t+y) = 



Then 



Rq(v) 



OO 00 

— OO mOO 






( 72 ) 



( 73 ) 

V 

( 74 ) 



( 75 ) 
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where 



f(C^,C2;y) 



- (^1 ^ 

^ Ci-Px") 

■ — 5 ^r~v~ 

cj 2tt (1-P,^)^ 



X 



X 



(76) 



and 



^X " 



Rx(y) 

R^CO) 



(77) 



Noting that the dependence of R upon y is only via p„j 

q A 

the notation R (P„) will be used for brevity. Then (75) 
q A 

becomes 



Rq(Px) 



N N 

i=i j=i ^ J 



X 



X 



X 



/ ^ / >5f(?^,C2iPx^^^l^^2 



i-1 ^J-1 



(78) 



Let 



01 - ?3_/ax 



(79) 



02 - ?2^*^X 



( 80 ) 



then ( 78 ) can be written as 



Rq(Px^ 



N N 
S S q q 
L=1 j=l J X^_3^ 



% i -(6,2-2pxPiB/e/) 

"x ‘’x — :;rT2^ 



/ / 

ibl 



"X '"X 



2(l-f^-) 



2(1 -p/)^ 



d0^d02 (81) 
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hIcvj 



but 



$1^ - 2B]_B2Px ^2 " ^^2 ■ 



Therefore 



- [ 



Bl^ - 2 PxPi^2 ■*■ ^2 

2(1 - p/) 



] = 



- (^2 - 

2(1 - Pjj2) 2 



and (8l)' becomes 



N N 



R (p„) =22 / “ 2 TT 



^ 2 
1 ^1/2 



Xj -(32 -Px3i)' 



q -X 



/ 



® 2 ( 1 - p ^ 2 ) 



i=l j=l X. 



1-1 



^ [2tt(1-p/)]'^ 
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Now 'let 



^2 - pA 

[2(1-Px2)3!s 



= a 



It follows that the second integral can be written as 



(erf[ - Px^l ]| 

^ [2(1-PX >r . [2(1-Px )] ' 



where 



2 

erf(a) = — |— / e dv 

ir^ 0 



( 82 ) 



(83) 



;l (84) 



(85) 



( 86 ) 



(87) 
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If a new function Is defined as 



N Qj <■ 
'p(Px»^l) ” ^ 2 I 



- Px^i 




] I (88) 



[2(l-p/)]^ 



then (8^) can be written as 




e 2 d6^ (89) 



2 

Dividing both sides of (89) by and substituting the 



expressions for X^ and q^, it is found that a new parameter 
appears, namely cr^/L . This parameter together with N and 
p^ will be the independent variables of (89). Note that with 
the introduction of d^/L , the dynamic range of the quantizer, 
L, is no longer an independent factor. 

By letting 



= q^/L 
X* = X^/L 



(90) 



(91) 



and 

it follows that 



(92) 



V 




rr * 



2 



and 
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♦'(Px.6i> = 





[2(1-Pjj2)]-5 




It can be shown, following a similar procedure, that for 

I 

p„ equal to one, (93) becomes 



autocorrelation of q. Equations (93) through (95) allow the 



a quantizer of N steps and saturation level L given that the 
Input is Gaussian, zero mean, and of known autocorrelation 
function. The procedure would be as follows. First, 
calculate j and q^ for the given quantizer and Input 



process. Second, find a value of for a particular y 
(Equation (76)). Third, form the function i|i'(p^,B^) for 
the particular value of p^^ determined in the second step. 
Fourth, replace the function i/^*(p^,B^) Into (93) and perform 
the integration. This procedure will map one point of the 
input autocorrelation function Into one point of the output 
autocorrelation function. 

A program was written to solve (93) through (95) for the 
case of sy^etrical quantizers (See Appendix under Transfer 
Characteristics of Odd level Quantizer, and Transfer Charac- 
teristics of Even level Quantizer) and the results are 




This gives us a 



2 



numerically so we can get the normalized 



q 



computation of the autocorrelation function of the output of 
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The 



plotted in Figures 8 through 1^1 as a function of Py. 

Rn(v) 

ordinate axis Is ^ ^ . In other words, for a given 

Oq2 q 

Py , P^ is obtained. The advantage of plotting it this way 
A q 

is that the curves are independent of the input autocorrela- 
tion function (as long as the input process is Gaussian and 
zero mean). 

In general, the curves are a function of (except for 

the special case of N=2). Note that for N=2 the result is 
the well-known relationship [Ref. 7] 



^q^^X^ 2 

— 2 — = ? Px 



(96) 



It is of interest to note the behavior of these curves as 
Oy/L approaches extreme values. For a quantizer with odd 
number of levels, as gets very large, the quantizer will 

appear as a hard limiter, hence the curves will approach that 
of Figure 8. As cfy/L approaches zero the quantizer will 
appear as a three step quantizer (eventually it will become 
an open circuit) and the curves will show a behavior similar 
to those shown in Figure 9. For a quantizer with an even 
number of levels, the behavior of the curves as becomes 

large is also as des.cribed above. When cTy/L approaches zero, 
however, the quantizer will appear again as a hard limiter, 
therefore the curves will approach that of Figure 8. 
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B. QUANTIZER FOLLOWED BY A SQUARE LAW DETECTOR 



The reason for analyzing this case will become apparent 
in Chapter V. Refer to Figure 6 and observe that . 



Y(t) = q^(t) 

and 

q(t) = Q[X(t)] 

Ryfy) = E[Y(t)Y(t+y)] 



(97) 

(98) 

(99) 



Upon substituting (97) and (98) into (99) it follows 



that 



Ry(p) = E{Q^[X(t)]Q^[X(t+y)]} 



( 100 ) 



Following a similar procedure as done for the previous case, 
(100) becomes 



Ry(Px) 






-3 



2 (^) f A(p^,3^) e 



(2tt)^ 1=1 ''^X 



dB- 



( 101 ) 



and 



-i N q! p I" 

A(py>3,) = ^ Z (^f|erf[ 



] - erf[ 



j=l X 



Xj_i/aj^ - Py3 

[2(1-P^2);,^ 



■'x^ 

-X 






( 102 ) 



For Pj{-=1 it can be shown that 



E[Y^] _ 1 ? 



I 2 (^)^ [erf ( ^ ) - erf (.• — ^, ) ] 

^ 1=1 ^X a’ (2)^ a' (2)'^ 






0.(2)^ 



(103) 
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defined in (90), (91), (92) and (77) 

respectively . 

A computer program was written for this case when the 

quantizer is symmetrical (See Appendix under "Transfer 

Characteristics of Odd Quantizer Followed by a Square Law 

Detector" and "Transfer Characteristics of Even Quantizer 

Followed by a Square Law Detector"). The same comments as 

in the previous section apply for the plots of this case 

Ry(vi) 

except that now the ordinate axis is p - - / - a - s - Instead of p„ 

u ; 1 

since Y(t) is not zero mean. The results are shown on 
Figure 15 for N equal to five and eleven. 

Ry(y) 

It can be shown that for the case of no quantizer ^ ^ y 
is given by 



Ry(p) 1 + 2py.^ 

RT^oy " 3 



(104) 



Note how the curves of Figure 15 depart from the curve given 
by (104) due to the effects of the quantizer. 

In summary, in this chapter an algorithm has been devel- 
oped-suitable for computer programming in order to find the 
autocorrelation function of the output of a quantizer (or 
quantizer and a square law detector) given that the Input 
process Is Gaussian and zero mean, and with known auto- 
correlation function. What It basically does (for both cases) 
is to reduce a double integral (75) and (100) Into a single 
integral by using the error function which is available In 
most computer libraries. 
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The transfer curves for both cases, shown on Figures 8 
through and 15 are Independent of the Input process 
providing It Is Gaussian and has zero mean. These results 
will be used In Chapter V. 








K-SCRLF.-4.00E-01 UNITS INCH. 
'f'-SCRLE-^. OOE-01 UNITS INCH. 



FIGURE 15 



^9 



( 



V. DEGRADATION FACTOR CALCULATION WITH 
BOTH SAMPLING AND QUANTIZATION 



In this' chapter the effects of sampling and quantization 
on the performance figure of the two digital radiometers 
under study are investigated. As in Chapter III, the two 
cases will be considered separately. 

A. IF SAMPLING 

1. Exact Results 

Refer to Figure 2. In Chapter III an expression was 

developed for the degradation factor for the case where there 

was no quantizer — see (30). Those results cannot be applied 

directly to this case since in order to obtain (30) it was 

assumed that q (kA) and q (kA) were samples of Gaussian 
c s 

processes — see (22). With the addition of a quantizer it 
is clear that this assumption no longer holds. Hence we 
must start again with a more general approach. This can be 
accomplished by concentrating first on the top channel of 
the radiometer of Figure 2. The processing that the ADC and 
computer do to the signal can be drawn sequentially as in \ 

\ 

Figure l6. It is argued that the order in v/hlch the signal 
q (t) is sampled and squared is immaterial as far as the 
random variable I^ is concerned. Hence Figure l6 is a valid 
model of a channel of the radiometer of Figure 2. 

Referring to Figure l6 it follows that [Ref. 12] 

^ (l-^)[Ry(kA) - E^CY]] (105) 

Xi n k=l n X 
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RyCk ) is the k— sample of the autocorrelation function of 

Y(t) 

and 



E[I^] = E[Y(t)] = Oq 



( 106 ) 



Kc(t) 




1 "Z 
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Using ( 26 ) and the fact that for the case of Figure 
2 the RF, Mixer, IF filter is symmetrical about f^, then 



4 = 2 al 



(107) 



and 



E[I] = 2E[I, ] = 2 



( 108 ) 



Applying the definition of degradation factor it 

follows that 

,2 



2a: 



p2 = 



„ n-1 , . 

“TP ip ^ — p ^ [Ry(kA) - a ] 

^ k=l T Y q 



da^ 

^c 



(109) 



“X ^*’’op 
c ^ 
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Using the same argument as in Chapter III (A<<t) it 



can be shown that the term 







V [Ry(ki) - 0^ ] 

k=l ^ ^ 



( 110 ) 



can be disregarded due to its small value compared with the 
other terms in the numerator of (110). 

Then (109) becomes 



= 



An 



Ac 



n-1 

E 

k=l 



[RY(kA) 



a** ] 

Q 



dcj 

^c 



at T 



op 



( 111 ) 



Looking at (111) term by term, it is noted that in 

the niimerator, the first term can be obtained from (103) and 

(95). The second term can be obtained from the corresponding 
RX^(kA) 

value of 5 77 rr~- and the transfer characteristics of a 

%c^°^ 

quantizer followed by a square law detector (Figure 15). 

2 4 

Actually, a scale factor of c„ /o„ will appear because of 

■^c 

the way the curves were normalized. The denominator can be 

obtained from (95). Note however that Equations (93), (95) 

eind the transfer characteristics developed in Chapter IV have 

Oy/L as a parameter and not cr„ . This is due to the choice 

^c 

of Cj^/L as a parameter since it was found in that chapter 
that the ratio of and L was of Importance and not their 
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absolute magnitudes. Therefore it makes sense to have (111) 
be a function of /L at T Instead of /L can 

A_ op A_ 

c ^ c c 

be thought of as a normalized standard deviation of X^(t) and 

its value can be related to temperature in the same form as 

Equations (10) and (11). 

2 2 2 

A plot of E[I]/L versus /L is shown in 

■^c 

Figure 17 for different values of N. The effects of satura- 
tion of the quantizer can be observed as temperature (equiva- 
2 2 

Idnt to a„ / L ) is Increased. 

■^c 

The four filters previously considered for this case 
were used again and the results are shown in Figures l8 
through 4l. The bottom curve for each figure represents the 
case investigated in Chapter III where there is no quantizer. 
The programs used to evaluate (111) are Included in the 
Appendix. It should be pointed out that the programs used 
to calculate the transfer characteristics of Chapter IV did 
not provide a continuous function as It Is shown on the plots 
of that chapter. Since In order to solve (111) the transfer 
curves must be known at every point, a least -mean-squared- 
error- polynomial fit was done to the data obtained from the 
programs of Chapter IV. 

Looking at Figures l8 through 4l It Is noted that In 

each case there Is a reasonable optimum value of /L 

•^c 

where the degradation factor Is lowest. Physically this can 

be explained by the fact that for small /L , the step 

•^c 

size degrades the performance of the system, hence Increasing 

Oy /L will tend to minimize this effect. However a point Is 
■'^c 
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reached where saturation effects start dominating and the 

further Increase of /L causes an Increase of the 

degradation factor. For N > 4 , /L optimum is around 

■^c 

0.35 regardless of the sampling rate. This would imply that 

for a given T , the gain setting of the RP, Mixer IF filter 

should be set at a value such that /L is equal to 0.35 

■^c 

for optimum performance. Note that this optimality is for 

small source temperatures. For large signals, dynamic range 

\ 

is a problem. See Chapter VI. 

2 . An Approximation 

The evaluation of (111) for the exact value of the 
degradation factor requires extensive computer time. For a 
quantizer of five bits or greater the computation time per 
case is very long (on the order of hours). Therefore it is 
of interest to investigate an approximate solution for the 
degradation factor and to compare the results with the exact 
values already obtained. 

Figure 42 shows one of the channels of the radiometer 
of Figure 2. Assume that the sampling rate is not too fast, 
the step size of the quantizer is small with respect to the 
standard deviation of the process X (t) and no saturation 

r 

takes place in the quantizer. Then 

q^(kA) = X^(kA) + (112) 

t h 

where X (kA) is the k— sample of X (t) and e, is a uniformly 

C C xC 

distributed random variable that accounts for the error 
introduced by the quantizer. The probability density 
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function of is the following [Ref. 13] 








Then 



( 113 ) 



1 

n 



n 2 
2 q^(kA) 
k=l ° 



1 

n 



n 

E 

k=l 



[X^OcA) + 2ej^X^(kA) + e^] 



(114) 



and 



E[I^] 



1 

n 



I E[X^(kA)] 
k=l ° 




Z E[2e^X^(kA)] +- Z E[e^] 
k=l k=l 



(115) 




\ 



FIGURE 42 



80 



Assuming e, and X(kA) are uncorrelated, which Is a 

kv C 

good assumption If the step size Is small compared with a . 

■^c 



Edi] ^ + 




( 116 ) 



Now 



E[I^] = -^ E [ S X^(jA) 
n j=l 



+ 2 S e.X^(jA) 
J=1 ^ ° 



+ 



n 



E 

J=1 



] • 



[ E xF(kA) + 2 E e, X^(kA) E ej ] 

kFl ° k=l ^ ° k=l ^ 



(117) 



Expanding (117) and using (23) and the assumption 

that e, and e. are Independent for k j , which Is good If 
K J 

the sampling rate Is slow. It follows that 





n 



H ^ 2 la 

i E (n-k)Rx (ki) + i y 

n k=l c 




H 

a 

720 



(118) 



It was shown previously that for the case of Figure 
2, If H^(f) was symmetrical about f^, then 



and 




(119) 



E[l] = 2 ECI^^] 



( 120 ) 



Using (119) and (120) It follows that 



2B^ n-1 2 J 1 

= BgA + -ij— I (1 - -) (kA) + ^ BgA + -ip ^ V . (121) 

Ov k=i c Ov 



But the first two terms are just the degradation 
factor for the case of no quantizer, therefore (121) can be 
written as 



= ^NQ ^ ^ 



B„A + 



a 



B„A 



3^0 



( 122 ) 



where Pj^q is the degradation factor with no quantizer. The 
last two terms of (122) represent a correction, due to 
quantization, for the degradation factor with no quantizer. 

The approximate solution is very good if the assump- 
tions under which it was derived hold. For N=17, /L = 0.3 

■^c 

(a/a^ = 0.42) and WA= 1.0 the error is less than one percent, 

■^c 

For values of a /L greater than 0.3 the approximation is not 

■^c 

good due to saturation. For values much smaller, it again 

falls, due to the large value of a/o^ . Figure 4 3 shows how 

■^c 

the two results compare for N=7 and different values of 

/L. In it, the behavior of the two 'solutions as /L is 
c .. : c 

increased can be observed. For N much greater than 7» the 

approximation is very good when /L is near 0.3* 
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B. RC FILTERING THEN SAMPLING 



Refer to Figure 3* Unfortunately there Is no exact 
solution for this case since there Is no way to get the 
second order probability density function of the process 
z(t). This would be needed to approach the problem In the 
same way as done for the radiometer of Figure 2. Two 
approximations will be discussed. 

1. First Approximation 

It was shown In Chapter III that 



and 



E[Z] = g(N^/2) 



4 = 



^Vrc 



(123) 



(124) 



Owing to the Central Limit Theorem It can be assumed 
that z(t) Is a Gaussian process If >> 1 . Due to the 

square law detector. It Is obvious that z(t) can never be 
negative, hence the assumption Is good only If the mean of 
z(t) Is sufficiently large compared to Its standard 
deviation so that the assumed Gaussian probability density 
function Is ■ negligible for negative arguments. 

It Is also assumed that the quantizer has a bias ’b* 
as shown on Figure 44. When b=0 this quantizer reduces to 
the one shown In Figure ?• It Is argued that the order In 
which z(t) Is sampled and quantized Is Immaterial as far as 
the output I Is concerned. Hence Figure 45 Is a valid model 
of the processing done by the ADC and computer. 
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FIGURE 44 
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FIGURE 45 



Let 

E[Z] = m (125) 

then the probability density function of z(t) is 

(z-m)^ 

f (z) = e (126) 

^ (2-nalr 

It can be shown, using a procedure similar to that 
Introduced in Chapter IV that 

E[q(t)] 






’4(p2j3) 



N q. 

Z ^ [ erf(— ^ 
1=1 ^ (2)% 



X. T-m 

-) - erf(^^) ] 



Z 

Xj^-m 



(2)%, 



N 



-B^/2 



(2tt)'^ i=l ^ 



Z q^L / iKP 2 »B) e dg 






(127) 



(128) 



X^-m 

N q^ C a„ " 

erf[ ^ 



= . if 

J=1 ^ 



P.B 



[2(l-p2)]'^ 



- ] - erf [ 



- p 6 

^7^ Pz^ 



[2(1-P^)]^ 



] 



} 



(129) 
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But from Figure 



= 2L(i - |)/(N-1) + b 



( 130 ) 



and 



= 2L[i - (N+l)/2]/(N-l) 



(131) 



Equation (130) can be written as 




(132) 



where X is the value of X. when the quantizer is zero 

Oi i 

centered. 



through (129), the parameter b-m appears. Physically, this 
is the offset of the mean of the process z(t) with respect 
to the center of the quantizer. When no source is present 
the value of this offset will be determined by the quiescent 
conditions of the radiometer. When looking at a source, 
however, the value of the offset will change since the mean 
of z(t) changes linearly with temperature. 



Substituting (131) and (132) into Equations (127) 



Let 



f = m 



o 



b 



X133) 



where m^ is the expected value of z(t) at Then (127) 

through (129) can be written as 
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SC2m=4 , ^(erfC 

^ ^ 1=1 oX '• 



N q» . - f’ +(2Bt^c^-^ (az^ - oz) 



(2) ^a* 



(134) 



- erf[ 



^1-1 - 



(2pa 






] 



V^’ 



Rq(Pz) _ 



N 

E 



^Z’ 






(2ir)'^ 1=1 ^Z 



Ci-f 



i|>'(p2,3) e 



-e^/2 



d3 



(135) 



(P2»P^ 



E I 

J=1 2a* t 



o' ■ PZ^ 
erf [ ^ p ] 

[2.(l-<yp]^ 



x<!,i.i-^* 



- P2^ 



- erf [ 



[2(1^^)] 



-'} 



(136) 



where o„ Is the value of the standard deviation of z(t) at 
T^p, and all primed variables Imply that they have been 
normalized by L. 

For the case of p^ = 1 (135) becomes 




N 

E [erf ( 

1=1 202 



^1 - 
a^(2) ■ 



- erf ( 





)] 



(137) 



Equations (134) through (137) can be solved In the 
same form as done In Chapter III. However, the addition of 
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two new parameters, (f* and (2Btj^Q)'^) , makes the problem 
difficult to solve in a general way due to the number of 
variables involved. 

One way to get around this problem is to assume that 
(2Bgtp^)'^ is large. Then 



N 

E[q] = Z q. 

i=l ^ 




f^CzjmjO) dz 



where 



f^CzjmjO) 



1 

(2Tra^)^ 



(z-m)^ 



2a‘ 



(138) 



(139) 



From (123) and (12^0 
m = 



Therefore 

N df2(z,m,a) dz 

dm ” . *^1 dm 

1--'- 

but 

^ 3f + ^ 

dm ” 3m 3a dm 



(140) 



(141) 



V 



(142) 



Substituting into (l4l) and using the fact that 



!£z 

3m 



if 

3z 



(143) 
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It follows that 






N 

Z q^[f2(X^_^,m,a)-f2(X^,ni,a)] 



N 

E 

1=1 






X. 



3a 



da (144) 



1-1 



' t 

If (2Bgtj^^) Is large. It means that the expected 
value of z(t) Is much greater In value than Its standard 
deviation. Therefore It can be assumed that a change In 
temperature affects the expected value of z(t) but the 
effects on a^ are negligible. In that case the second term 
of (l44) can be disregarded. 

Using (66) It follows that 



!l‘= + 2 a 



na, 



n-1 

E 

k=l 



(1 - 



R (kA) - E^[q] 
^)[ 2 ^ 



(145) 



Substituting (l44) and (145) Into (13), It Is easily 
shown that 






p2 = 



Or 



2 t 



RC 



. n-1 R (kA) - E^[q] 

2B^t E [ -S ^ ] 

^^^^RC k=l 



( Mal)2 

^ daz ^ 



2Bg(A/tp^)^ n-1 R^(kA) - E^[q] 



T/t. 



RC 



E k [ -3- 
k=l 



f dE[q] ,.2 
^ daz ^ 



(146) 
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The parameter 2Bgtp^^ on the numerator of (l46) will 
cancel out. The third term of (l46) can be neglected only 
if (see comments in Chapter III. The addition of a 

quantizer does not invalidate them) . 

u 

It- must be pointed out that the parameter (2Bgtpj^)^, 
even though it does not appear in (l46), is an important one. 
It is a measure of how much smoothing the RC filter does. 

A large value makes quantization difficult. On the other 
hand, a small value invalidates the assumption that z(t) is 
a Gaussian process. Other consequences of this parameter 
will be discussed in Chapter VI. 

Equation (l46) can be evaluated using (134) through 
( 137 ) and (l44). Basically it is done in the same form as 
for the case discussed in Section A of this chapter. Using 
( 135 ) and ( 136 ) a set of transfer characteristics has to be 
computed similar to Figures 8 through l4 (that set of curves 
is the solution when f'=0). The major difficulty lies in 
the additional parameter f’ which increases the dimensionality 
of the problem. 

The case when f'=0 (no offset) was analyzed using 
Figures 8 through l4 and (l44). The results are shown in 
Figures 46 through 51. For any other value of f’ (134) 
through ( 137 ) have to be evaluated. 

As an example the case of N=13, C2/L=0.2 and f’=0.5 
was solved. The results are shown in Figures 52 and 53* 

Note that they are very close to the zero offset case for a 
seven level quantizer with a^^/L^O . The reason is that 
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a/a 2 is the same for both cases. They are not exactly 
equal however, since in one case clipping occurs essentially 
for positive values of the signal, while in the other case 
clipping occurs symmetrically for both positive and negative 
values of signal. Figure 52 is a plot of the transfer char- 
acteristics for this example. The programs used are in the 
appendix (See "Transfer Characteristics of Quantizer With 

Offset," and "RC Filtering Then Sampling, Offset Case"). 

\ 

The case when N=2, f*=0 is of interest since results 
are obtained easily without extensive computer programming. 
Using (1^6) and ( 96 ) it can be shown that 



p2 _ IT A . A 

^ “ IT t — ^ t — 
^ ^RC ^RC 



n-1 

Z 

k=l 



arcsin 



R2(kA) 



^RC 




n-1 

E 

k=l 



k arcsin 



R2(kA) 



(147) 



The last term can be neglected when T>>tj^^. As the 
sampling rate Increases, A goes to zero and Equation (147) 
becomes 



F^ 



/ arcsin e~^ dy 
0 



(148) 



For large t ,-.the upper limit can be replaced by 
infinity and (l48) integrates to 

F^ = In 2 (149) 

therefore 

F = 1.04 ( 150 ) 
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The results are surprising since it means that a 
hard limiter degrades the system only by four percent. 
However, the dynamic range of a radiometer using a hard 
limiter would be very small hence this result is of little 
practical use. More on the problem of dynamic range will 
be discussed in Chapter VI. 

2. Second Approximation 

Using the same arguments as in part A of this 
chapter, it can be shown that [Ref. 9] 



P 



2 




a 



0 = 21 . 



(151) 



The results obtained from (151) are in close agree- 
ment with those obtained using the first approximation even 
for low values of N. It must be kept in mind that the 
second approximation does not take into account saturation 
effects of- the quantizer. For the case when f’=0 it was 
found that if cr^/L < .6 WA > 1 , ^ N > 5, 

the second approximation was within 1.5 percent of the first 
(For'N=17 and the above conditions, the difference was less 
than 0.6 percent). These results indicate that a sizable 
amount of saturation can occur before (151) no longer holds. 
Figure 53 shows how the two approximation compare for the 
example done on page 91. Figure 5^ shows a comparison for 

N=7. f*=0 and different values of a /L. The different 

Z 

behavior of the two approximations as Increased can 

be observed. 
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To summarize. In this chapter an exact solution has 
been found for degradation factor of the radiometer of 
Figure 2. For the radiometer of Figure 3 no such solution 
exists, but an approximation can be made under the reasonably 
good assumption that z(t) Is a Gaussian random process. For 
both cases the major difficulty In global analysis lies In 
the number of parameters Involved. 

An approximate solution was developed which Is In 

\ 

close agreement with the first solution. Its major advantage 
Is that It Is easy to compute and as the number of steps 
Increases, It gets better. For practical ADC It should be 
adequate. 
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VI. DYNAMIC RANGE AND LINEARITY OF A DIGITAL RADIOMETER 



In this chapter the effects of the addition of an ADC 
on the linearity of theoutput of a radiometer is investigated. 
The main problem here is the saturation effect of the 
quantizer. 

A. IF SAMPLING 

Refer to Figure 2. The output voltage versus temperature 
characteristics for this case was investigated previously, 
and is shown in Figure 17. If the ADC was ideal (no satura- 
tion and a/(J = 0) there would be one straight-line curve 
with slope equal to two. It can be seen that for small 

/l 2 (proportional to temperature) the curves tend to have 
-^c 

a slope of two. 

Figure 17 provides all the information needed for the 
problem under investigation in this chapter. However, the 
way the curves were normalized makes them awkward to use. 

They can be changed starting with (108) and (95). 




) - erf ( -, - ^'V ] (152) 



x: 



but 




c 



2 



L 



(N-l)a/a 



(153) 



X 



c 
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Therefore (152) becomes 



E[I] _ 






N 

E 

1 = 1 ' 



q ’ (N-1) 2 , x’ (N-1) 
Co^-ZT-] jerfC ^ 



■2oJ7K' 






2(2)%^/a 



] - erfC 



X^l^CN-l) 
2(2) V /a 



] (154) 



Equation (154) is plotted in Figure 55- The curves asymptote 
to 



E[I] _ (N-1)' 

„2 " 2 



(155) 



This can be shown by using the fact thas as 
gets large the quantizer appears as a hard limiter. 

Figure 55 shows that a substantial range of linear 

» 

output can be obtained with this scheme before saturation 
dominates . 

In order to use Figure 55 a criterion must be 

defined for how much departure from linearity is allowed 

for the temperature range expected (one percent, for example). 

Once that is defined. Figure 55 can be entered to determine 

the number of steps of the quantizer required. The value of 

the abscissa at T will be determined mainly by how much 

op 

degradation of performance is allowed (Chapter V). We thus 
have a trade-off of linear dynamic range versus minimum 
detectable temperature AT. 
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B. RC FILTERING THEN SAMPLING 

Refer to Figure 3- A similar analysis can be made for 

this case. However the additional parameters (f* and 

h 

(2BEtRc) ) make a general representation difficult. 

I 

Basically, (127) has to be rewritten in terms of cr^/a 
instead of It can be shown easily that (127) becomes 



N q’ r 
E[I] = a" 1 {erf[ 

a 1=1 2sal! ( 






(2)"^ s al 



] - 



erf[ 



X^-1 - f + (q^^-<rg)s(2B^t^g) 
(2)*^^s a" 



( 156 ) 



where f ’ , and are as previously defined, and 






(157) 



P ff _ 






(158) 



s = 2/(N-l) 



Figure 56 Is a plot of (156) for the case of f’=0, 

u 

=0.5 and (2BT:.tT3n) = 1.0. The abscissa axis Is shown 

Lq £j nO 

as a function of '^'^'^op^* curves asymptote 

^o ^ 

to (N-l)/2. Some observations are of Interest. First, the 

h 

slope near the origin Is equal to (2Bgtj^^) . Since for 

a given quantizer E[I]/a reaches a fixed value for large 
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values of , increasing (which physically 



o 



means that more smoothing will be done by the RC filter) 

decreases the range of linear output. Therefore t^^ .should 

be chosen as the minimum value which will enable the sampler 

to operate properly. Second, the parameter will be 

^o 

determined by the amount of degradation tolerated at . 

And third,. the desired range of linear output will determine 

the number of steps of the quantizer. Therefore a procedure 

to use a Figure such as 56 would be as follows: First, 

determine and cr” using the criteria given above. 

^o 

Second, draw a set of curves of the expected value of the 
output versus temperature using (I56). Third, decide upon 
a criterion of how much departure from linearity is allowed 
in the temperature range of operation, and finally obtain 
the required number of steps from the curves drawn. 

To summarize, the effects of the addition of an ADC on 
the linearity of the output of the radiometers under 
consideration, have been investigated in this chapter. The 
knowledge of the behavior of the output is important for 
calibration purposes. 
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9 



values of cr^/cr^ , Increasing (2Bp,tpp) (which physically 



o 



E"RC' 



means that more smoothing will be done by the RC filter) 

decreases the range of linear output. Therefore tj^j^ .should 

be chosen as the minimum value which will enable the sampler 

to operate properly. Second, the parameter a!^ will be 

^o 

determined by the amount of degradation tolerated at 

And third,. the desired range of linear output will determine 

the number of steps of the quantizer. Therefore a procedure 

to use a Figure such as 56 would be as follows: First, 

determine (2Bgtj^j^) and o” using the criteria given above. 

^o 

Second, draw a set of curves of the expected value of the 
output versus temperature using (I 56 ). Third, decide upon 
a criterion of how much departure from linearity is allowed 
in the temperature range of operation, and finally obtain 
the required number of steps from the curves drawn. 

To summarize, the effects of the addition of an ADC on 
the linearity of the output of the radiometers under 
consideration, have been investigated in this chapter. The 
knowledge of the behavior of the output is Important for 
calibration purposes. 
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where P is the degradation factor for the equivalent total 
power radiometer. The factor multiplying P is the performance 
figure of a balanced Dicke radiometer. 




PIGURE 57 



Equation (l60) follows from the fact that for the half 
of the switching cycle that the radiometer is connected to 
the source^ it behaves as a total power radiometer. 

Care must be taken however on the selection of x. If 
commencement of sampling is delayed after switching, it must 
be corrected for the loss in integration time. Also, for 
the case of RC filtering followed by sampling, the introduc- 
tion of the RC filter must be reflected in an Increase of x 
since P takes only into account the effects of the inclusion 
of the ADC. 

The problem of dynamic range and linearity of the output 
can be handled with the curves developed in Chapter VI. 
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B. NOISE ADDING RADIOMETER 

Figures 59 and 60 show two digital noise adding radio- 
meters. The one In Figure 59 corresponds to the IF sampling 
case and only one channel Is shown. The one in Figure 60 
corresponds to the RC filtering case. The extension of 
results already found to NAR Is considerably more difficult 
than for a balanced Dlcke radiometer. This Is due to the 
processing Involved [Ref. 13]* 

Let 



I = E + e 
n n n 



( 161 ) 



and 



I = E^ + e 
o o o 



( 162 ) 



where I and I^ are the outputs of the radiometer with 
no 

T _+T and T ^ respectively. E^ and E^ are the expected 
op n op n o 

value of the output for the case of T +T and T respec- 

^ op n op ^ 

tlvely. and are independent, zero mean random variables 

with standard deviation equal to a_ and Oj , and represent 

■^n "^o 

the deviation from the mean for a particular measurement 
(one cycle). 



Then the Y-factor is [Ref. l4] 
I E 1 + ^ 

Y = j2. = g2. [ ] 

0 0 Gq 



(163) 
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If e„ and are small compared with E and E , then 
no ^ no 



K e 

Y = — ^ [ 1 + — ^ 



n 



( 164 ) 



E[Y] = 



n 



(165) 



and 



Yon o 



( 166 ) 



Using Equation (l66) and the fact that for a NAR 



= (1 + _2£) ( JL. 

op n E 



(167) 



where t is the switching period; for the case of Figure 2 
it can be shown by straight substitution that. 



I c ( 



lECln] 2 

d4 ’ 



n‘ 



NAR 



X 2 dE[I ] 
(1+-52-) f 2 + ( ^ 



I'd ^ 



:°E- 



(1 + 



,2 

r 



n 



ECln] 



dE[I^] 



da; 



T +T 
n op 



ETin 



n- 



•op 



(168) 



where a„ is a„ at T , a„ is a„ at T 

X„ X op ' x„ X n . 

o n 
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FIGURE 60 



is the degradation factor of the corresponding total 

power radiometer at '^”'^n'*’'^op * ^o degradation 

factor of the corresponding total power radiometer at T=T^p. 

2 

For the case of Figure 3, must be replaced by in (168). 
All the terms of (l68) can be calculated from equations 
previously developed. 

It must be noted that for the case when the signal is 

filtered before sampling (Figure 60), f (offset) cannot be 

zero at T ^ or T +T since the ADC would give zero output, 
op op n 

The problem of dynamic range and linearity is easy to 
handle. From Figures 55 or 56, depending on the system under 
consideration, and (165) a plot can be made relating E[I] 
to temperature. Figure 6l shows how this is done for the 
system of Figure 59* 

In summary, in this chapter the results obtained in the 
six previous chapters for the digital total power radiometers 
shown in Figures 2 and 3, have been used to determine 
performance of other digital radiometers. The balanced 
Dicke and noise adding radiometers were investigated. For 
the balanced Dicke radiometer, there is a simple and straight- 
forward relationship. For the NAR, however, a relationship 
was found but it is more involved. The problem of dynamic 
range and linearity can be handled easily. 



117 



STC STO ZTO CQS 9SD COO 000 




K-SCflLE=3.00E+00 UNITS INCH- 
V-SCRLE-3.00E+00 UNITS INCH. 



FIGURE 6l 



118 



VIII. DIGITAL FILTERING 



In the analysis of the classical total power radiometer 
an Integrator is used to smooth the output signal. This 
scheme works well but it has the disadvantage that it gives 
output information only at discrete Intervals of time, corre- 
sponding to blocks of data. Such performance might not be 
acceptable in some cases because a smooth output might be 
desirable. Therefore, a pure Integrator is not generally 
used but rather a low pass filter (such as a RC) that pro- 
vides a continuous output. The theory and equations of a 
classical total power radiometer do not change when using 
these low pass filters instead of an Integrator if an 
equivalent integration time is defined [Ref. 1]. 

The digital total power radiometers of Figures 2 and 3 
suffer the same disadvantage as the classical TPR. That is, 
output information is obtained only after the summers have 
completed the summation of the n samples. Therefore a scheme 
where output information is obtained more often must be 
investigated. This chapter deals with the substitution of 
the summers of Figures 2 and 3 by digital filters. There is 
an additional reason for wanting to do this analysis. Since 
a computer is already being used in the radiometers of 
Figures 2 and 3, the use of a digital filter often does not 
require any more hardware. 

Figure 62 shows a digital filter operating on a digital 
signal Xj^. The filter is characterized by its Impulse 
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FIGURE 62 

response and the output is The filter is assumed 

to be causal (i.e. hj^=0 for k<0) and time invariant. This 
filter can replace the summers of the radiometers defined in 
the previous chapters. 

Referring then to Figure 62 



Y 



k 



00 



I X h 

i=_oo 



k-i 



(169) 



By straightforward algebra it can be shown that 



E[Y. ] = E[Y] .= E[X] E h 
^ i=0 



and 



= E h. E [Ry - E^[X]3h^ 
^ 1=0 ^ m=0 m-i 



( 170 ) 



(171) 



Where R^^ = E[Xj^Xj^^^] 



is the discrete autocorrelation of Xj^. 



Let 

j = 1 - m 



( 172 ) 
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then 



Oy = Z 

^ i=0 



1 p 

h. Z [Ry - E^[X]]h, . 

^ j=-«> 



(173) 



Assume that h decays much slower than • This can be 

J 

seen by observing the location where this filter would go 
in Figures 2 and 3* For the case of Figure 2, the power 
spectral density at the input of the digital filters is 
shaped by the RF, Mixer, IF filter which is broad-band. 

For the case of Figure 3> it can be assumed only when tj^^^ 
is small (this is the case for which Figures 46 through 51 
apply ) . 

Then (173) becomes 



" Jo j=-«» (174) 

Let be the value where H(jQ) is negligible. Then 
(174) can be rewritten as 

^ o i “ i 

aj = Z hj- Z iRy - E^tX^] + Z h^ Z [R^ - E^[X]] 

^ i=0 ^ j=-co i=j^ 1 j=-oo 

(175) 

It can be shown that the first summation is much smaller 

than the second, and hence it can be ignored. Furthermore 

the limits of the second summation can be changed as follows 

2 

without altering the value of Oy significantly 
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OO 00 

Oy = Z E [R^ - E^[X]] (176) 

1=0 J=“~ j 



For the summer used In previous chapters. It can be shown 
(using the same notation of Figure 62) that 



22 ^-^ 2 k ? 

Oy ^ Oy + 2n E [R^ - E^[X]]- 2n E ^ [R^ - E^[X]] (177) 

^ ^ k=l k=l ^ \ 



and 

E[y] = E[X] (178) 

The upper limit of the summation of (177) can be replaced 
by Infinity (due to the narrowness of Ry,(^)) and the third 
term can be Ignored (see Chapters III and V). 

Substituting (170) and (176) Into (4) for the general 
digital filter. It follows that 



op 



®X ” 2' 

E hj + 

1=0 

op 



00 00 

E h^ E [Ry 



E^[X]] 



( Eh / (^)2 

i=0 ^ 



(179) 



op 



Doing the same with (177) and (178) It follows that 



op 



22 _ 



^ "^op^ k=l 






E^CX]] 



r dE[X]- ,2 
L jjip J 



( 180 ) 



T 

op 
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Now if a comparison is made of (179) and (l80) it is 
clear that they are equal if n^ is defined as follows 




00 



( ^ h 

i=0 



Z 

i=0 



h 



2 

i 



( 181 ) 



and ng can be interpreted as an equivalent summation number. 
Physically it is a measure of the performance of a general 
digital filter in terms of the performance of a pure summer. 
In other words, it provides a way to replace the digital 
filter by a summer of n^, samples without changing the per- 
formance of the radiometer. Then for a given sampling rate 
and known n^ an equivalent integration time can be determined 
and a performance figure can be calculated using the curves 
of Chapter V. 

For a digital total power radiometer such as those of 

AT 

Figures 2 or 3 using a digital filter, ^ — would be given by 

op 

^ ^ r P (182) 

•^op (Bj.ngA)'^ 

where F is the degradation factor for the equivalent 
radiometer using summers and (B^n^A) ^ is the performance 

figure of a total power radiometer with Integration time 
equal to n^A . 
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It Is often convenient to express (l8l) in its equivalent 
z-transform since the transfer function of a digital filter 
is usually given in terms of its z-transform. It can' be 
shown that [Ref. 15] 



(H(z) 












H(z)H(z”^)z”^dz 



(183) 



where the contour integration denoted by f is around a unit 
circle and H(z) is the transfer function of the filter. 

■ As an example of the use of (I 83 ) a first and second 
order digital filter will be considered. 

For a first order filter, the transfer function is 
given by [Ref. 13] 



H(z) = : ^ 

1 - rz 



It can be easily shown that (I 83 ) becomes 




1 H- r 

1 - r 



(184) 



(185) 



The transfer fundtlon of a second order filter is given 
by 

H(z) = ^ 5 ( 186 ) 

1 - z” a - z” X 

Substituting (I 86 ) into (I 83 ) and after considerable algebra, 
it can be shown that 
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(187) 



- X^(a^+2) - 2Xa^ - + 1 

" (l-a-X)2 (,_,2) 



Equation (182) and (I 83 ) are very similar to the ones 
defining the equivalent integration time of a low pass filter 
of an analog total power radiometer. This is not surprising 
since many similarities exist between digital and analog 
networks even though they are not equivalent . 
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IX. CONCLUSION 



This work investigates the performance of radiometers 
using digital processing. It is assumed that the changing 

I 

of the signal from an analog to digital form is performed 
before the final output is obtained. Two digital total 
power radiometers were considered in detail. 

In Chapter II a figure of merit is defined in order to 
evaluate performance. The effects of the addition of an 
analog to digital converter are discussed. 

In Chapter III the effects of an ideal ADC are investi- 
gated. It serves as an introduction to the problem of 
sampling and quantization. The results are useful since 
they provide a lower bound for the degradation factor as 
defined in Chapter II. 

Chapter IV deals exclusively with the effects of a 
quantizer on the autocorrelation function of a zero mean 
Gaussian input process. The problem of a quantizer followed 
by a square law detector is investigated. 

Chapter V combines the information of Chapters II, III 
and IV and treats the effects of a real ADC. Both quantiza- 
tion and sampling effects are investigated. For one of the 
radiometers under investigation, exact results are obtained. 
For the other a reasonably good assumption is made and an 
approximate solution is obtained. Other approximation are 
discussed and compared. 
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Chapter VI discusses the problem of dynamic range and 
linearity involved in a digital radiometer. 

In Chapter VII the results obtained in the previous 
chapters are used to investigate other types of digital 
radiometers. Balanced Dicke and noise adding radiometers 
are considered. 

Chapter VIII investigates the problem of using a general 
digital filter for processing of signals in a radiometer. 

The major difficulty encountered throughout the work was 
the number of variables that enter into the problem. There- 
fore it was impossible to do a global analysis. However, 
the factors and considerations that enter into the design of 
digital radiometers were treated in sufficient depth so that 
a good understanding of the advantages, problems and 
limitations of digital radiometers is presented. 
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APPENDIX 



The following is a list of the programs used to 
evaluate the equations and draw the figures throughout the 
thesis. They are listed in the order they were used. 



TRANSFER CHARACTERISTICS OF ODD LEVEL QUANTIZER 
N=NUMBER OF STEPS 

VAR=NORMALIZED STANDARD DEVIATION OF X 
RHC=COPRELATION COEFICIENT OF X 
DEL=INTEGRATION STEP 
XF=X(N ) 

XC=X(0) 

IMPLICIT REALMS! A-H,C-Z) 

DIMENSION SYS(5001),0(34),X(34) , FUN (5001) ,W(18,2) 

1,M(18.2) 

C PROBLEM PARAMETERS 

READ(5.2C7)N 
207 F0RMAT(I2) 

WRITE{6.666)N 
666 FORMAT (• ',12.//) 

DO 205 NN=1,8 
READ(5.200)VAR 
200 FORMATiFA.l) 

C ALGORITHM PARAMETERS 

CEL=.001 
RH0=0.0 
R00=0.0 
XF=1000. 

X0=-1000. 

SY=-5.0 
SS=0. 0 
S=0.0 
. E=0.0 

STP=(2./(N-1.) ) /VAR V 

PP=.3989A23 
M = (N-1 )/2 
L=DABS( SY) /DEL 

C QUANTIZER CHARACTERISTICS 

DC 81 1=1, N 

0( I ) =( I-(N+1.1 /2. 0)--^STP 
X( I ) = ( I-N/2.0)=PSTP 
81 CONTINUE 

C LIMITS OF INTEGRATION 

CO 97 J = 1,L 
SYS( J) =SY 
SY=SY+DEL 
97 CONTINUE 

DO 91 IR=1,N1 
IF(IR.NE.l) GO TO 92 
M(l,l)=l 

K3 = ( X( l)+5. 0) /DEL+1. 

W0=X(1)-SYS(K3) 

M( 1,2) = K3 
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W(1 ,2 )=WG 
GO TO 91 

92 N(IR,1)=M(IR-1,2)+1 
W(IR.l)=DEL-WO 
K4=(X(IR) + 5.0) /DEL+1. 

WO=X( IR )-SYS (K4 ) 

IR,2)=K4 
W( IR ,2 )=WO 
91 CONTINUE 
C FINAL POINT 

DO 98 1 = 1, N 
IF(I.NE.l) GO TO 99 
E2=X0/l.-f.lA 
E1=X( I )/ l.-^14 
GO TO 100 

99 IF(I.NE.N) GO TO 101 
E2=E1 

E1=XF/1.414 
GO TO 100 

101 E2=E1 

E1=X( I ) /1.414 
100 R2t9FRF( E1)-DERF(E2) 

RR=0( I )=^=«2=i'R2 
SS=‘SS+RR 
98 CONTINUE 

RATT=SS/2.0 
RAT=DSORT( RATT) 

WRITE (6, 102)RATT,"RAT 

102 F0RMAT(2F15.9) 

SY=-5.0 

DO 30 MH=1,21 
C FORM FUNCTION PHI 

A = 2.0*( 1.0-RH0=^=^2) 

AA=DSORT ( A) 

DO 2 J=1,L 
DO 1 1=1, N 
IF(I.NE.l) GO TO 5 
E2=(X0-RH0*SY) /AA 
E1=(X( I )-RHO«SY)/AA 
GO TO 25 

5 IF( I.NE.N) GO TO 3 
E2=E1 

E1=(XF-RH0^SY)/AA 
GO TO 25 
3 E2=E1 

E1=(X(I )-RHO^SY)/AA 
25 R1=DERF( E1)-DERF(E2) 

R=OI I )=^=R1 
S=S+R 

1 CONTINUE 

FUN( J) = { S/2.0 )*OEXP( I SY=«=>}=2)/ (-2.01 ) 

. S=.0 
SY=SY+DEL 

2 CONTINUE 

C MAIN INTEGRATION 

DO 7 I1=1,N1 
44 IF( II .NE.l ) GO TO 11 
M1=M( 1,1) 

M2=M(1,2) 

CORR=W( 1 ,2)=^^FUN(M2 ) 

GO TO 15 
11 M1=M(I1,1) 

M2=M( II ,2) 

CORR = W( I 1, l)=f=FUNIMl) + W( II ,21*FUN(-M2> 
15 M3=M2-1 

DO 21 K0=M1,M3 
011=FUN(KG)+FUN(K0+1) 

B=B+011 
21 CONTINUE 

R0=0( I 1)^( DEL*B+CORR) 

FCO=ROO+RO 

B=.0 
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7 CONTINUE 
POOO=ROC^PP 
POOO=(1.0/RATT)*ROOC 
WRITE (6»40)RHD ♦ROOD, ROQO 
40 FORMAT (3 D2 5. 9) 

WRITE! 7»204)RQ00»R0CC 

204 FORMAT! 2 F8. 5) 

RCQ=0.0 

S=.0 

SY=-5.0 

IF!MM.E0.20) GO TO 67 
RHO=RHO+.05 
GO TO 30 
67 RH0=RH0+.03 
30 CONTINUE 

205 CONTINUE 
STOP 

END 



TRANSFER CHARACTERISTICS OF EVEN LEVEL QUANTIZER 
N=NUMBER OF STEPS 

VAR=NORMALIZED STANDARD DEVIATION OF X 
RHG=C0RRELATI0N COEFICIENT OF X 
DEL=INTEGRATION STEP 
XF=X!N) 

XO=X 10) 

IMPLICIT REAL=^-8!A-H»0-Z) 

DIMENSION SYS! 5001) »0!34) ,X!34) t FUN! 5001) ,W!18,2) 
1 ,M!18»2 ) 

C PROBLEM PARAMETERS 

N=4 

WRITE!6,666)N 
666 FORMAT! • • » I2»//) 

DO 205 NN=1»8 
READ!5,200) VAR 
200 FORMAT !F4.1) 

C ALGORITHM PARAMETERS 

RHO=0.0 
CEL=0.001 
R0Q=0.0 
XF=1000. 

XC=-1000. 

SY=-5.0 

SS=0.0 

S=0.0 

B=0.0 

STP=!2./ !N-1.) 1/VAR 
PP-. 3989423 
Nl=N/2 

• L=CABS!SY)/DEL+1 
C' QUANTIZER CHARACTERISTICS 

CO 81 1 = 1, N 

0! I) = ! I-lN + 1 .)/2.0 )=^STP 
Xm=! I-N72.0)=rSTP 
81 CONTINUE 

C LIMITS OF INTEGRATICN 

DO 97 J=1.L 
SYS! J)=SY 
SY=SY+DEL 
97 CONTINUE 

DO 91 IR=1,N1 

IF !I R.NE.l ) GO TO 92 

M! 1, 1) =1 

K3=!X!l)+5.0)/DEL+l. 

WO=X! 1 )-SYS! K3 ) 

M! 1,2)=K3 
W!1 ,2 ) =WC 
GO TO 91 

92 M! IR, 1 )=M! IR-1 ,2)+l 
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W( IR ♦! )=DEL-WO 
K4=(X( IR) + 5.0) /DEL + 1 . 

WO=X( IR )-SYS(K4) 

M( IR,2 ) =K4 
W( IR ,2) =W0 
91 CONTINUE 
C FINAL POINT 

DO 98 1=1, N 
IF( I.NE.l) GO TO 99 
E2=XC/1.414 
E1=X( I ) /I. 414 
GO TO 100 ' 

99 IF(I.NE.K) GO TO 101 
E2=E1 

E1=XF/1.414 
GO TO 100 

101 E2=E1 

E1=X( I )/1.414 
100 R2=DERF( E1)-DERF(E2) 

RR=0(I )-*^:^2i=R2 
SS=SS+PR 
98 CONTINUE 

RATT=SS/2.0 
RAT=DS0RT(RATT ) 

V^RITE(6,102)RATT,RAT 

102 F0R,MAT(2F15.9) 

SY=-5.0 

DO 30 MM=1,21 
C FORM FUNCTION PHI 

A=2.0i'(1.0-RH0**2) 

AA=DSORT(A) 

CO 2 J=1,L 
DO 1 I=1,N 
I F( I.NE.l) GO TO 5 
E2 = (X0-RH0-e-SY)/AA 
El=( X( I )-RHO*SY) /AA 
GO TO 25 

5 IFU .NE.N) go to 3 
E2=E1 

E1=(XF-RH0*SY) /AA 
GO TO 25 
3 E 2 Ei 1 

ll = {X{ I )-RHO=^SY)/AA 
25 R1=DERF(E1)-DERF(E2 ) 

R=0( I )AR1 
S=S+R 

1 CONTINUE 

FUN{ J) = ( S/2.0)=^-DEXP( ( SY=«‘*2 ) / (-2 . 0 ) ) 
S=.0 

SY=SY+DEL 

2 CONTINUE 

C - MAIN INTEGRATION 
DO 7 I1=1,N1 
44 IF( I I.NE.l) GO TO 11 
M1=M(1.1) 

M2=M( 1.2) 

CORR=Wl 1,2)*FUN(M2) 

GO TO 15 
11 M1=M{ 11,1) 

M2=M{ 11,2) 

C0RR=W(I1,1) 4:FUN(M1 ) + W( 1 1 , 2 )*FUN ( M2 ) 
15 M3=M2-1 

CO 21 KQ=M1,M3 
011=FUN{K0)+FUN( KO+1 ) 

8 = 8+011 
21 CONTINUE 

R0=0{ I 1)^( DELi'B+CORR) 

RO0=ROO+F:O 

E =.0 

7 CONTINUE 
ROOO=RGO*PP 
R000=( 1 .0/RATT )^ROOC 
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WRITE(6»49) RHOtROeCtROQO 
40 FORMAT (3D25. 9) 

WRITE(7,204) ROOOtRQCO 

204 FGRMAT( 2F8. 5) 

R00=0.0 

S=.0 

SY=-5.0 

IF(MM.E0.20) GO TO 67 
RHC=RHC+.05 
GO TO 30 
67 RH0=RHO+.03 
30 CONTINUE , 

205 CONTINUE 
STOP 
END 



TRANSFER CHARACTERISTICS OF ODD LEVEL QUANTIZER 
FOLLOWED BY A SQUARE LAW DETECTOR 
N=NUMBER OF STEPS 

VAR=NORMALIZED STANDARD DEVIATION OF X 
RHC=CORRELATIGN COEFICIENT OF X 
DEL=INTEGRATION STEP 
XF=X(N) 

XC-=X ( 0 ) 

IMPLICIT REAL*8( A-H,C-Z) 

DIMENSION SYS( 5001 ) tO( 34) , X( 34) ,FUN ( 500 1 ) , W ( 1 8 ,2 ) 
1 tM(18,2 ) 

C PROBLEM PARAMETERS 

N=17 

CO 205 NN=1,2 
READ( 5,200) VAR 
200 FORMAT! F4.1) 

WRITE (6 ,888)VAR 
888 FORMAT! • • ,F4. 1) 

C ALGORITHM PARAMETERS 

RHC=0.0 
DEL=.001 
R0Q=0.0 
XF=1000. 

xn=-iooo. 

SY=-5.0 

SS=0.0 

S=0.0 

B=0.0 

STP=!2./!N-1.) ) /VAR 
PP=. 3989423 . 

Nl=!N-l)/2 
L=DABS! SY) /DEL 

C QUANTIZER CHARACTERISTICS 

■ DO 81 1=1, N 
Q!I)=!I-!N+1.) /2.0)*STP 
X! I ) = ! I-N/2.0)*STP 
81 CONTINUE 

C LIMITS OF INTEGRATION 

CO 97 J=1,L 
SYS! J) =SY 
SY=SY+DEL 
97 CONTINUE 

DO 91 IR=1,N1 
IF!IR.NE.l) GO TO 92 
M!l,l )=1 

K3 = !X! D+5.0) /DEL+1. 

WO=X! 1)-SYS!K3) 

M! 1 ,2) =K3 
W! 1,2)=W0 
GO TO 91 

92 M!IR,1)=M!TR-1 ,2)+l 
W! IR,1)=DEL-W0 
K4=!X! IR)+5.0 )/ DEL+1 . 
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WO=X( tP )-SYS(K4> 

M( IR,2 )=K4 
W( IP »2 )=lifO 
91 CONTINUE 
C FINAL POINT 

DO 98 I=1,N 
IF( I.NF. 1) GO TO 99 
E2=X0/1.A1A 
E1=X(I )/1.414 
GO TO 100 

99 IF(I.NE.N) GO TO lOl 
E2=E1 

E1=XF/1.414 
GO TO 100 

101 £2=E1 

E1=X( I ) /1.A14 
100 R2 = DERF( El )-D=RF(E2 ) 

RR=0( I )**4«R2 
SS=SS+RR 
98 CONTINUE 

RATT=SS/2. 0 
RAT=DSORT(RATT ) 

WRITE(6»102)RATT,RAT 

102 FORMAT! 2F15. 9) 

SY=-5.0 

DO 30 MM =1.21 
C FORM FUNCTION LAMNA 

A=2.0*(l .0-RH0^=?2) 

AA=DS0RT(A) 

DO 2 J=1.L 
DO 1 I=1.N 
IF(I.NE.i) GO TO 5 
E2=(XC-RH0=^=SY) /AA 
E1 = (X( I )-RHOT-SY)/AA 
GO TO 25 

5 IF( I.NE.N) GO TO 3 
E2=E1 

El = ( XF-RHO^SY) /AA 
GO TO 25 
3 E 2 — E 1 

El=( X( I 1-RHO^SY) /AA 
25 R1=DERF( El 1-D5RF(E21 
R=0( I )*=<'2*R1 
S=S+R 

1 CONTINUE 

FUN! J)=( S/2.01*DEXP( (SY**2>/ (-2.0) ) 
S=.0 

SY=SY+DEL 

2 CONTINUE 

C MAIN INTEGRATION 

CO 7 11=1, N1 
44. IF( II .NE.l ) GO TO 11 
M1=M(1,11 
M2=M(1,2) 

CORR=W( 1 ,2)*FUN( M2) 

GO TO 15 
11 M1=M(I1,1) 

M2=M( 11,2) 

CORR=W{ 1 1, 1)=^FUN(M1 )+W( 1 1,2)*FUN( M2) 
15 M3=M2-1 

DO 21 K0=M1,M3 
011=FUN(K0)+FUN (KO+1 ) 

B=B+011 
21 CONTINUE 

RQ=0( II )*^2*(DEL^B+C0RR) 

ROO=FOO+RO 

B=.0 

7 CONTIMUE 
ROOO=ROQv'PP 
R000=( 1. C/RATT)-ROOO 
WRITE(6 ,40 ) RHO, RCOC.POOO 
40 FORMAT ( 3D25. 9) 
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WRITE(7,230)R0Q0 
230 FCRMAT(F8.5) 

R00=0.0 

S=.0 

SY=-5.0 

IF(MM.E0.20) GO TO 67 
RHC=RH0+.05 
GO TO 30 
67 RH0=RH0+.03 
30 CONTINUE 
205 CONTINUE 
STOP 
END 



TRANSFER CHARACTERISTICS OF EVEN LEVEL QUANTIZER 
FOLLOWED BY A SQUARE LAW DETECTOR 
N=NUNBER OF STEPS 

VAR=NORMALIZEO STANDARD DEVIATION OF X 
RHO=CORRELATICN COEFICIENT OF X 
DEL=INTEGRATION STEP 
XF=X(N) 

X0=X(0) 

IMPLICIT REAL^i'BI A-H,0-Z) 

DIMENSION SYS( 5001 ) ,0( 34 ) » X ( 34 ) ,FUN ( 500 1 ) , W ( 1 8 ,2 ) 
1 ,M(18»2 ) 

C PROBLEM PARAMETERS 

N=4 

DO 205 NN=1»3 
RE AD (5, 2 00) VAR 
200 F0RMAT(F4.1) 

WRITEI6 ,888 )VAR 
888 FORMAT! • • ,F4. 1) 

C ALGORITHM PARAMETERS 

RH0=0.0 
DEL=.001 
P0Q=0.0 
XF=1000. 

X0=-1000. 

SY=-5.0 

SS=0.0 

S=0.0 

B=0.0 

STP=(2./(N-1.) )/VAR 
PP=. 3989423 
Nl=N/2 

L=CA8S(SY)/0EL+1 

C QUANTIZER CHARACTERISTICS 

DO 81 1 = 1, N 

.. 0(I)=(I-(N+1.) /2.0)-^STP 
X( I )=( I-N/Z.OI ^-STP 
81 CONTINUE 

C LIMITS OF INTEGRATION 

DO 97 J = 1,L 
SYS( J)=SY 
SY=SY+DEL 
97 CONTINUE 

DO 91 IR=1,N1 

IF( IR.NE.l ) GO TO 92 

M( 1, 1) = 1 

K3=(X(1 )+5.0)/0EL+l. 

WO=X( 1)-SYS(K3) 

M( 1,2)=K3 
W(1,2)=WC 
GO TO 91 

92 M( IR, 1 )=M( IR-1 , 2 )+l 
W(IR,l)=DEL-WO 
K4=( X{ IR )+5. 0) /OEL+1. 

WO=X( IR)-SYS(!<4) 

M( IR,2 ) =K4 
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W( IR ,2)=W0 

91 CONTINUE 

DO 91 IR=1,N1 
IF(IR.NE.l) GO TO 92 
M ( 1 , 1 )=1 
DO 93 K3=l ,L 
W0=SYS(K3)-X( 1 ) 

IFIWO.GT .0 .0 ) GO TO 9A 

93 CONTINUE 

94 M(1,2)=K3-1 
W(i,2)=DEL-W0 
GO TO 91 

92 M( IR,i )=M( IR-1 ,2)+l 
W(IR,1)=W0 

DO 95 K4=1,L 
WO=SYS (K4)-X(IR ) 

IF(WO.GT.O.O) GO TO 96 

95 CONTINUE 

96 M(IR,2)=K4-1 
W( IR ,2)=DEL-W0 

91 CONTINUE 
C FINAL POINT 

DO 98 I=1,N 
IF(I.NE.l) GO TO 99 
E2=X0/1.414 
E1 = X( I ) /I. A14 
GO TO 100 

99 IF(I.NE.N) GO TO 101 
E2=E1 

E1=XP/1.414 
GO TO 100 

101 E2=E1 

E1=X( I )/1.414 
100 R2=DERF( E1)-DERF{E2) 

RR=0( I 
SS=SS+RR 
98 CONTINUE 

RATT=SS/2.0 
RAT=DSORT( RATT ) 

WRITE! 6, 102)RATT,R>AT 

102 FORMAT (2F15. 9) 

SY=-5.0 

DO 30 MM =1,21 
C FORM. FUNCT ION LAMNA 

A=2.0*-(1.0-RH0^v2) 

AA=DSORT(A) 

CO 2 J=1,L 
DO 1 1=1, N 
IF(I.NE.l) GO TO 5 
E2=(XO-RHO^SY)/AA 
E1=I X( I )-RHO=^SY) /AA 
GO TO 25 

5 IF(I.NE.N) GO TO 3 
E2=E1 

E1=(XF-RH0*SY) /AA 
GO TO 25 
3 E2=E1 

E1={X{ I )-RHO*SY)/AA 
25 R1=DERF( El )-DERF( E2 ) 

R=Q( I )*:*2?R1 
S = S+R 

1 CONTINUE 

FUN( J) = ( S/2.0)*DEXPn SY^*2 ) / ( -2- 0 ) ) 
S=.0 

SY=SY+DEL 

2 CONTINUE 

C MAIN INTEGRATION 

DO 7 11=1, N1 
44 IF(Il.NE.l) GO TO 11 
M1=M(1 ,1 ) 

M2=M( 1,2) 

CORR=W( 1,2)^FUN(M2) 
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GO TO 15 
11 M1=M(I1,1) 

M2=M( 11,2) 

CORR=W( II ,1 )tFlJN(Ml )+W ( 1 1 , 2 ) ’*«FUN ( M2 ) 
15 M3=M2-1 

DO 21 K0=M1,M3 
011=FUN( KO)+FUN( KO+1) 

B=B+011 
21 CONTINUE 

R0=0( I1)=^:*2’^(DEL*B+C0RR) 

ROQ=ROO+RO 

B=.0 

7 CONTINUE 
ROOO=RCO=?^PP 
RCOC=( 1 .0/RATT )*RQQO 
WRITE! 6,40) RHO.RQGC ,RGQO 
40 FORMAT! 3025. 9) 

WRITE!?, 230)R00C 
230 F0RMAT!F8.5) 

R00=0.0 
S = .0 
SY=-5.0 

IF!MM.E0.20) GO TO 67 
RHO=RHC+.05 
GO TO 30 
67 RHC=RH0+.03 
30 CONTINUE 
205 CONTINUE 
STOP 
END 



C EXPECTED VALUE OF OUTPUT AND DERIVATIVE I.F. SAMPLING 

C N=NUMBER OF STEPS 

C XF=X!N) 

C X0=X!0) 

IMPLICIT REALMS! A-H,0-Z) 

DIMENSION Q!34) ,X!34) , Y!440) ,SIG!440) ,V!440) ,Z(440) 

XF=100.0 

XC=-100. 

READ!5,71)N 
71 FORMAT! 12) 

S=.0- 
VARR=.02 
CO 8 L=1 ,440 
VAR=DSORT! VARR) 

STP=!2./!N-1.) )/VAR 

CO 1 1=1, N 

X! I)=! I-N/2. )=?STP 

0!I) = ! I-!N + 1.)/2.)=^STP 

IF!I.NE.l) GO TO 5 

E2=X0/1.414 

E1=XI I )/1.414 

GO TO 6 

5 IF!I.NE.N) GO TO 3 
E2=E1 

E1=XF/1.414 
GO TO 6 
3 E2=E1 

E1=X!I ) /1.414 

6 R1=DERF!E1)-DERF!E2) 

R=Q! I )*^2*R1 

S=S+R 

1 CONTINUE 
OARR=VARR*S 
Y!L)=OARR 
V!L)=VARR 
SIG!L)=S/2.0 
VARR=VARR+.01 
S=0.0 
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8 CONTINUE 

CALL DDET5( .01 , YtZ ,AAO, lER ) 

WRITE(6,82)(V(L),Y(L),Z(L),SIG(L) ,L=1,440) 
82 FORMAT! • • ,4F15.9) 

STOP 

END 



C DEGRADATICN FACTOR GF IDEAL FILTER 

IMPLICIT REALMS! A-H»0-Z) 

REALV8 ITITLE(12)/' S WETT • » 1 • •/ 

REAL’^'4 X3,F5»F6,F7»F8,F9.F10 
REALMS TITLE! 10)/10>^' '/ 

REAL^}:4 LABEL/4H / 

DIMENSION Z!43 ) .F2 !43) tWI !43) »Y !43-) , DELY ! 43 I , B ! 2 1 ) , 
1F!9.31),F5!31),F6!31),F7!31),F8!31),F9!31),F10!31), 
2 API! 31) ,AP2!31 ) , AP3 !3 1 ) » AP4 ! 31 ) , AP 5 ! 31 ) , AP 6 ! 31 ) , 
3F12!31) ,F13! 31) »SB!21) »X3!31) 

DO 98 1=1,22 

98 READ!5,99) WI !I ) 

99 F0RMAT!F7.1) 

DO 55 MM=1,6 

READ!5,999) VARl , VAR2 , VAR3 , VAR4 , V AR5 
999 FORMAT! 5 F4.1) 

CO 5 M=2,6 
DO 1 1=1,22 
READ!5,2)F2!I ) 

2 F0RMATIF8.5) 

1 CONTINUE 

READ!5,47)AA,F,C,D 

47 FORMAT !4F15. 7) 

A=AA-C^*2 

WRITE!6,48)A,E,C,D 

48 FORMAT! • • ,4F15.7) 

Z!21)=.98 

Z!22 ) = 1 .0 
S=0.0 

DO 3 1=1,20 
Z! I)=S 
S=S+.05 

3 CONTINUE 

C FIT PCLYNOMIAL 

CALL LS0PL2!22,-16,Z,F2,WI , Y ,DELY , B , SB , T ITLE ) 

E!1)=0.0 

WT=100. 0 

XX=.01 

DO 35 1=1,31- 
N=WT/XX-1. 0 
RUN=0.0 
DO 4 J=1 ,N 

IFU.GT.650) GO TO 49 
T = 3.1416’:^J^XX 
X7=DSIN!T)/T 
X=DABS! X7) 

IF!X.LT.0.1) GO TO 3332 

PUN=B! 1 )+B!2)*X+3!3 )=^X^=^2+B !4 )* X>f^»3 + B ! 5 ) *X*^4 

1- fB!6 )^X*=7=5 + B!7)=^X! 6 )+B ! 3 )*X=i--^7+B ! 9 ) =f^^X ! 8 ) + B ! 10)*X**9 

2- fB ! 1 1 )=i'-X=i=i‘10+3 ! 12 ) *X^- 1 1+B ! 13 )=^X=!^--*1 2 + B ! 14 ) =i<X''5=*13 

3- fB! IB)’' X*^14-»-3! 16)vX**15+B!17)*X*=^'16 
GO TO 3333 

3332 PUN=!F2 !3)-F2!l ) )*X’^10.0 

3333 RUN=RUN-hPUN 

4 CONTINUE 

49 A1=A*XX 
A2=2.0’^XX=«=E*RUN 
F4=2 .0=^! A1-hA2)/D=^^*2 
F!M,I)=DS0RT!F4) 

IF!I.NE.l) GO TO 66 
XX=XX-h.09 

GO TO 35 
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66 XX=XX+.l 
35 CCNTINUE 
5 CONTINUE 
XX=.l 

F(1 «1 )=1 .0 
X3( 1)=0. 0 

C NO QUANTIZER CASE 

DC 7 1=2.31 
N=WT/XX-1. 0 
R=0.0 

DO 8 J=1,N 
T=3.1416=^J*XX 
S=(DSIN(T 1/T)*=^2 
R=R+S 

8 CONTINUE 

YY=XX+2 .0*XX=^-R 
F( l.I )=DSORT(YY) 

X3( I )=XX 
XX=XX+.l 
7 CONTINUE 
WRITE(6.987) 

987 FORMAT! '1') 

WRITE(6.31 ) (X3 ( I ) . F(l. I ) .F(2, I) , F(3, I) , F(4, I ) ,F( 5, I ), 
1F(6.I).I=1,31) 

31 FORMAT (• ',7F12.7> 

WRITE (6. 888) 

888 FORMAT! ' * . /////) 

X=0.0 

READ!5.88)N 

88 F0RMAT!I2) 

WRITE!6. 89)N 

89 FORMAT! • • ,12) 

C APPROXIMATION 

DO 297 J=l,31 
RR = F! 1. J)5<‘^2 

SS=X*! 2.0/!N-l ) )*=^2/6. 0 , 

SSS=X^!1 .0/360 .0)-^! 2.0/ (N-li 
APP=RR+SS/VAR1^"*2+SSS/VAR1’^=?4 
API! J)=DSORT! APP) 

APP=RR+SS/VAR2=P=^=2 + SSS/VAR2=^^4 
AP2! J)=DSORT! APP) 

APP=RR+SS/VAR3’^-2 + SSS/VAR3**4 
AP3U) =DSQRT! APP) 

APP=RR + SS/VAR4«^2+SS5/VAR4=?^4 
AP4!J)=DS0RT!APP) 

APP=RR+SS/VAR5'^=^2+SSS/VAR5=^'^4 

AP5!J)=DS0RT!APP) 

X=X+. 1 

297 CONTINUE 

WRITE!6.il00)!X3!J).F!l,J),APl!J),AP2!J),AP3!J),AP4!J) 

1 ,AP5! J) , J = 1 .31 ). 

1100 FORMAT! • '.7F12.7) 

CO 18 K=l,31 V 

F5!K)=F!1,K) 

F6!K )=F! 2.K) 

F7!K)=F!3.K) 

F8!K)=F!4.K) 

F9!K)=F! 5.K) 

18 F10!K)=F!6.K) 

CALL DRAW! 31 . X3 . F5 . 1 . 0 , LABEL , IT IT LE , 0 . 5 , 1 . 0 , 0 , 0 , 2 , 2 , 6 , 

16. 1. LAST ) 

CALL DRAW!31 ,X3. F6,2.0,LABEL. ITITLE.0.5 , l.O.O.O, 2. 2, 6, 

16.1. LAST) 

CALL DRAW! 31.X3.F7.2. 0 . LAB EL , I T I TLE , 0 . 5 , 1 . C, 0.0, 2, 2, 6, 

16. 1. LAST) 

CALL DRAW! 31. X3.F3. 2.0. LABEL . I TITLE ,0. 5 ,1. 0,0, 0,2 ,2,6, 

16.1. LAST ) 

CALL DRAW! 31 ,X3 . F 9 . 2 , 0 . LABEL , I T IT L E , J . 5 , 1 . 0 . 0 , 0 , 2 , 2 , 6 , 

16. 1. LAST ) 

CALL DRAW!31 ,X3 . FIO . 3. 0 . L ABL 5. ITITLE,0.5,1.C,0,0,2,2,6 

1.6. 1. LAST) 

55 CONTINUE 
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STOP 

END 



C DEGRADATION FACTOR OF SINGLE TUNED FILTER 

IMPLICIT REAL*8( A-H,0-Z) 

REAL*4 X3»F5»F6tF7,F8tF9,F10 

REAL=^=8 ITITLE(12)/' SWETT',11*' */ 

REAL*8 TITLEd J)/10^' V 
REAL*4 LABEL/4H / 

DIMENSION Z(43)tF2(43) »WI(43),Y(43) ,DELY(43),B(21) , 
1F(9»31 ) tF5 (31) »F6(31 ) » F7 (31 ) tF8 (31 ) , F9(31) , F10I31 ), 
2APK31) »AP2(31) »AP3(31) tAP4(31) ,AP5(31) ,AP6(31) f 
3F12(31 ) » F13(31 ) »SB( 21) »X3( 31) 

DO 98 I=lt22 

98 READ(5»99)WI(I ) 

99 FORMAT (F7.1) 

DO 55 MM=lt6 

REA0(5»999 ) VAR 1, VAR 2»VAR3tVAR4, VARS 
999 FORMAT (5F4.1 ) 

DO 5 M=2t6 
CO 1 1=1,22 
READ(5,2)F2( I ) 

2 FORMAT(F8.5) 

1 CONTINUE 

READ! 5,47)AA,E,C,D 

47 F0RMAT(4F15.7) 

A=AA-C**2 

WRITE! 6,48)A,E,C,D 

48 FORMAT!' *,AF15.7) 

Z!21)=.98 

Z! 22) = 1.0 
S=0.0 

DO 3 1=1,20 
Z! I ) = S 
S=S+.05 

3 CONTINUE 

C FIT POLYNOMIAL 

CALL LSCPL2!22,-16, Z, F2 , W I , Y , DEL Y , B , SB , T I TL E ) 

B!1)=0.0 

WT=10.0 

XX=.01 

DO 35 1=1,31 
N=WT/XX-1.0 
RUN=0.0 
DO 4 J=1,N 

IF!J.GE.2C0), GO TO 4 
T=3. 1416=^XX=^ J 
X=DEXP!-T) 

IFIX.LT.0.1) GO TO 3332 

PUN=B!1 )+B!2)^X+P!3 )=='X’t:?-2 + B !4 )=^X^^3 + B (5 )*X**4 
1+0!6)^X**5+B! 7)-X^«^6+B ! 8)=^X=^*7+ B ( 9 ) X**8 + B ( 1 0 )=^^X=?^*=9 
2 + B!ll )^X^'^10+B! 12 )*X**1 1+S ! 13 ) ^X*=!= 12+B ( 14)’«'X**13 
3+B! 15)vX*^14+B ! 16)*X*-n5 + B(17)*X=^*16 
GO TO 3333 

3332 PUN= !F2 (3 )-F2! 1 ) )=?^X*10.0 

3333 RUN=RUN+PUN 

4 CONTINUE 
A1=A=?XX*3.1416 
A2=2.0»'3. 1416^ XX^E=^ RUN 
F4=2.0-^( A1 + A2) /D=1'^2 
F!M,I)=DSORT!F4) 

IF(I.NE.l) GO TO 66 
XX=XX+.09 

GG TO 35 
66 XX=XX+.l 
35 CONTINUE 

5 CONTINUE 
XX=.l 

F!1 ,1 )=1 .0 
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X3( 1) =3. 0 

C NO QUANTIZER CASE 

DO 7 1=2,31 
R=3. 1416*XX 
YY=R/DTANH( R) 

F(l,I)=DSORT(YY) 

X3(I)=XX 
XX=XX+. 1 
7 CCNTINUE 
WRITE! 6, S87) 

987 FORMAT ( • 1* ) 

WRITF( 6,31) (X3( I ) ,F(1 ,I ) ,F(2, I) ,F(3, I) , F(4, I ) , F(5, I ), 
1F(6, I ) ,1=1,31) 

31 FORMAT!' *,7F12.7) 

WRITE!6,888) 

888 FORMAT!' ',/////) 

X=0.0 

READ! 5, 88) N 

88 FORMAT! 12) 

WRITE!6,89)N 

89 FORMAT! ' ' ,12) 

C APPROXIMATION 

DO 297 J=l,31 

RR=F!1, J)=*^*2 

SS=X^=!2 .0/ !N-1 ) ) **2/6.0 

SSS=X* !1 .3/363.3)=*= ! 2.0 /!N-1) ) 

APP=RR + S S/ VAR l'**2+SSS/ VAR 1**4 
AP1!J)=DS0RT!APP) 

APP=RR+SS/VAR2**2+SSS/VAR2**4 
AP2! J)=DSORT!APP) 

APP=RR+SS/VAR3**2+SSS/VAR3*=^4 
AP3! J)=DSORT!APP) 

APP=RR+SS/VAR4**2+SSS/VAR4**4 

AP4! J) =DS0RT!APP1 

APP=RR + SS/VAR5*=i'2+SSS/VAR5**4 

AP5!J)=DS0RT!APP) 

X=X+.l 

297 CONTINUE 

WRITE!6,1100)!X3!J),F!1,J),AP1!J),AP2!J),AP3!J),AP4!J) 
1,AP5! J) , J=l,31 ) 

1100 FORMAT!' '.7F12.7) 

DO 18 K=l,31 
F5!K )=F! 1,K) 

F6!K) = F!2,K) 

F7!K)=F!3,K) 

F8!K)=F!4,K) 

F9!K) = F!5,K) 

18 F10!K) =F !6 ,K) 

CALL DRAW! 31, X3, F5, 1, 0, LABEL , I TITLE ,0. 5, 1.0, 0,0, 2, 2, 8, 

16.1 , LAST ) 

CALL DRAW! 31 ,X3, F6 , 2 , 0 , LABEL , IT IT L E , 0 .5 ,1.0, 0,0, 2, 2, 6, 

16. 1, LAST ) 

■ CALL DRAW!31,X3,F7,2,0,LABEL, ITITLE,0.5, 1.0,0,'3,2,2,6, 

16, 1, LAST) 

CALL DRAW! 31,X3, F8,2, C, LABEL , ITITLE ,0.5,1. 0,0, 0,2, 2, 6, 

16, 1, LAST) 

CALL DRAW!31,X3,F9,2,0,LABEL,ITITLE,0.5 ,1.0, 0,0, 2, 2, 6, 

16.1, LAST ) 

CALL DRAW! 31 ,X3 , FI 0 , 3 , 3 , L ABL5 , I T I T L E , 0 . 5 , 1 . 0 , 0 , 0 , 2 , ^ , 6 

1 .6. 1, LAST) 

55 CCNTINUE 

STOP 

END 



C DEGRADATION FACTOR OF GAUSSIAN FILTER 

IMPLICIT REALMS! A-H,G-Z1 
REAL*4 X3,F5,F6,F7,F8,F9,F10,F11 
REAL*8 T ITLE ! 10 )/10*» '/ 

REAL*4 LABEL/4H / 



DIMENSION Z(43),F2(A3) ,WI(43),Y(43) ,DELY(43) ,8(21) , 
1F(9,31),F5(31),F6(31),F7(31),F8(31),F9(31),F1D(31), 
2APK31) ,AP2(31 ) ,AP3(31) ,AP4(31) , APS ( 3 1 J , A P6 ( 3 1 ) , 
3F12(31),F13(31),SB(21),X3(31),F11(31) 

DO 98 1=1,22 

98 READ(5,99)WI(I ) 

99 FGRMAT(F7.1) 

DO 55 MM=1,2 

READ( 5, 999 )VAR1, VAR2,VAR3 
999 F0RMAT(3F4.1 ) 

DO 5 M=2,4 
DO 1 1=1,22 
READ(5 ,2 )F2 ( I ) 

2 FORMAT(F8.5) 

1 CONTINUE 

READ (5 ,47) AA,E,C,D 

47 FORMAT! 4F15.7) 

A=AA-C^-*2 

WRITE(6,48) A,E,C,D 

48 FORMAT! ' • ,4F15. 7) 

Z!21 )=.98 
Z!22)=1.0 

S=5.0 

DO 3 I =1 ,20 
Z! I )=S 
S=S+.05 

3 CONTINUE 

C FIT POLYNOMIAL 

CALL LS QPL 2! 22,-16, Z, F2 , WI , Y ,DELY , B , SB , T I TLE ) 

U=2.773 

B! 1)=0. 0 

WT=10.0 

XX=.01 

DO 35 1 = 1,30 
N=WT/XX-1.0 
RUN=3.0 
DO 4 J=1,N 

IF!J.GE.100) GO TO 4 
T=! !3.1416-XX«J)**2)/U 
X=DEXP!-T) 

IF!X.LT.0.1 ) GO TO 3332 

PUN=B! 1)+B!2)*'X+B!3 )?X»'^i2+B!4)^X*^3 + B! 5 ) ->=X’S-*4 

1 + B!6)^X-^^^5 + B!7 )*X^=^6+B! 8)~X«*7+B! 9)«X1=*8+B! 10)^X**9 

2 + R!ll )':■ X*^10+3 ! 12 )*X--^11+B ! 13 )’^'X**12 + B! 14 I^fX’f-IB 
3+B! 15)*X*^14+B! 1 6 )=» X^* 1 5+B ! 1 7 ) « X^* 16 

GO TO 3333 

3332 PUN=(F2!3)-F2! 1) )=)'X*10.0 

3333 RUN=RUN+PUN 

4 CONTINUE 
A1=A^XX*1. 505 
A2=2,0=^1.50 5TXX*E*RUN 

. F4=2 .0^(A1 + A2)/C*=!'2 
F( M,I ) =DS0RT!F4) 

IF(I.NF.l) GO TO 66 
XX=XX+.09 
GO TO 35 
66 XX=XX+.l 
35 CCNTINUE 

5 CONTINUE 
XX=.l 

F( 1,1) =1.0 

C NO OUANTIRZER CASE 

CO 7 1=2,30 
N = WT/XX-1. 0 
R=0.0 
X3!l)=0.0 
S=0. 0 

DO 8 J=1,N 

T=!2.0* !3. 1416*XX’!=J )**2 )/U 
S=DEXP!-T) 

R=R+S 

8 CCNTINUE 
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V= 1 . 505=^^ XX+ 3. 0 1* XX=«= R 
F(l, I) = DS0RT(V ) 

X3(I )=XX 
XX=XX+. 1 
7 CONTINUE 
WRITE (6«987) 

987 FORMAT! • 1' ) 

WRITE (6. 31 )(X3m,F(l»I),F(2,I),F(3,I),F(4,I),I = l,3 0) 
31 FORMAT! • • ,5F12.7) 

WRITE!6,888) 

888 FORMAT!' •«/////) 

X=0.0 

READ!5«88)N 

88 FORMAT! 12 ) 

WRITE!6«89)N 

89 FORMAT!' ',12) 

C APPROXIMATION 

DO 297 J=l,30 

RR=F!1, J )*-^2 

SS=X* ! 2 . 0/ ! N-1 ) ) ^=^=2 /6 . 0 

SSS=X*! 1.0/360.0)»! 2.0/!N-l) 

APP=RR + SS/VAR1T*2 + SSS/VAR1*>^4 
APKJ) =DSORT! APP) 

APP=RR + SS/VAR2-<*2+SSS/VAR2*^4 
AP2! J)=OSORT! APP) 

APP=RR+SS/VAR3=^=i=:2+SSS/VAR3*=P4 
AP3!J)=DS0RT! APP) 

X=X+.l 

297 CONTINUE 

WRITE!6,31)!X3!J),F!1,J),AP1!J),AP2!J),AP3!J),J=1,30) 
DO 18 K=l,30 
F5IK )=F! 1,K) 

19 F6!K) = F!2,K) 

20 F7!K)=F!3,K) 

21 F8!K) = F!4,K) 

F9!K)=AP1!K) 

F10!K)=AP2!K) 

18 F11!K)=AP3!K) 

CALL DRAW!30,X3 , F5 , 1 , 0, LABEL , IT IT L E , 0 . 5 , 0 . 5 , 0 , 0 , 2 , 2, 6, 

16. 1, LAST) ) 

CALL DRAW! 30 , X3 , F6 , 2 , 0 , LAB EL , I TI TLE , 0 . 5 ,1.0, 0,0, 2, 2, 6, 

16.1, LAST) 

CALL DRAW! 30 , X 3 , F7 , 2 , 0 , LABE L , I TI TLE, 0.5, 1.0, 0,0, 2, 2, 6, 

16.1 . LAST) 

CALL DRAW! 3 3 ,X3 , F8 , 2 , 0 , LAB EL , I T ITL E , 0 .5 ,1.0, 0,0, 2, 2, 6, 

16. 1. LAST ) 

CALL DR AW! 30, X3, F9, 2, 3, LAB EL ,ITITLE,0.5,1.C,0,0,2,2,6, 

16. 1. LAST) 

CALL DRAW! 30, X 3, FI 0,2, 4, LAB EL, I TITLE ,0. 5, 1.0,0, 0,2, 2, 6 

1.6. 1. LAST) ■ 

CALL DRAW! 30,X3 , FI 1 ,3 , 5 , LABEL , I TI TLE , 0. 5 ,1.0, 0,0, 2, 2, 6 

1,6, l.LAST) 

55 CONTINUE 
STOP 
END 



C DEGRADATION FACTOR C*" BUTTERWORTH FILTER 

IMPLICIT REAL=5=8! A-H.O-Z) 

REAL*4 X3, F5.F6, F7, F8,F9,F10 

REALMS ITITLEI12)/' SWETT',11*' '/ 

REAL^8 TITLE! 10) /lO*' '/ 

REAL*4 LABEL/4H / 

DIMENSION Z!43) , F2 ! 43 ) , W I ! 43 ) , Y !43 ) ,DELY!43) ,B!21 ) , 
1F!9,31),F5!31),F6!31),F7!31),F8!31),F9!31),F10!31), 
2AP1 !31 ) ,AP2 !31 ) , AP3 !31 ) , AP4!31 ) , AP 5 ! 3 1 ) , AP 6 ! 3 1 ) , 
3F12I31) ,F13!31) , SB!21) ,X3!31) 

CO 98 1=1,22 

98 READ!5,99)WI !I ) 

99 FORMAT! F7.1) 
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CO 55 KM=1,6 

READ(5«999) VARl , VAR2 » VA R3 , VARA , V AR5 
999 FORMAT ( 5 F4. 1) 

CO 5 M=2»6 
DO 1 I =1*22 
READ(5«2)F2(I ) 

2 FORMAT (F8. 5) 

1 CONTINUE 

RFAO( 5,47)AA,E,C«0 

47 F0RMAT(4F15.7) 

A=AA-C^'^ 2 

WRITE(6,48.)A,E»C,D 

48 FORMAT! • ' «4F15.7) 

Z( 21 )=.98 

Z(22 ) = 1 .0 
S=0.0 

DO 2 I=l«20 
Z( I )=S 
S=S+. 35 

3 CONTINUE 

C FIT POLYNOMIAL 

CALL LS0PL2(22,-16, Z.F2.WI »Y,DELY» B, SB t TITLE) 

B(1)=0.0 

WT=10. 0 

XX=.01 

DO 35 I=l«31 
N=WT/XX-1. 0 
RUN=0. 0 
DO 4 J=1,N 

iF(J.GT.30J) GO TO 49 
T=1.414'PDEXP(-J^3. 1416>'XX/1.414) 

X7=T*DC0S( J*3. 14 16*XX/ 1.414-3.1416/4.0) 

X=DABS( X7) 

IF(X.LT.O.l) GO TO 3332 

PUN=B( 1 )+B(2 )*X + e (3 )=^'X^=^=2+B(4)^X^*3 + B( 5 )*X"=«=4 
l+B(6)4^Xi'«5+S( 7)*X^*6+P( 8)*X*’«'7+B(9)TX^=48 + B (10 )«X=^*9 
2 + e(lU^X«=j=10+B( 12)=^X-^'^'11+B( 13)*X^=^-12+B( 1 4 ) =^'X'-i‘^ 13 
3+B(15)^X**14+B(16)=PX**15+B(17)-+X=i^*16 
GO TO 3333 

3332 PUN= (F2 (3)-F2( 1 ) )^X*10.0 

3333 RUN=RUN+PUN 

4 CONTINUE 

49 Al=A’f-'XX*1.48 

A 2 = 2 . 3 - 1 . 4 8 ^ - X X * E R U N 
F4=2.0'^( Al+A2>/D^*2 
F(M, n=DSQRT(F4) 

IF( I .NF. 1) GO TO 66 
XX=XX+.09 
GO TO 35 
66 XX=XX+.l 
35 CONTINUE 
5. CONTINUE 
XX=.l 

F(l.l)-1.0 

X3(1)=0.0 

U=1.414 

C NO QUANTIZER CASE 

DO 7 1=2,31 
N=WT/XX-1. 0 
R=0.0 

00 8 K=1,N 

0=3.1416*K!‘XX*U 

S = 0EXP(-0) t(1.0+CSIN(0) ) 

R = R+S 

8 CONTINUE 

YY = 1 .48i^XX+2 .9 6T-XX*R 
F( 1, I ) =DSORT(YY) 

X3( I )=XX 
XX=XX+.l 
7 CONTINUE 

WRITE(6,987) 

987 FORMAT! • 1' ) 
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WRITE(6,31)(X3(I),F(ltI)»F{2,I),F(3,I)>F(4,I),F(5,n, 
1F(6» I ) » I=l»31) 

31 FGPMATC *,7F12.7) 

WRITE( 6,888) 

888 FORMAT(» ',/////) 

X=0.0 

READ( 5,88) N 

88 FORMAT (I 2) 

WRITE(6,89)N 

89 FORMAK • » ,12) 

C APPROXIMATION 

DO 297 J=1.31 

RR=F( 1, J)=^=^2 

SS=X^(2 .0/ (N-1 ) )=«-*2/6.0 

SSS=X* ( 1 . 0/360 . 0 I'i' ( 2 .0/ ( N-1 ) ) ^*4 

APP = RR + SS/VAR1**2+SSS/ VAR 1^*4 

APK J)=0S0RT(APP) 

APP=RR+SS/VAR2*=^ 2+SSS/VAR2’^^4 
AP2( J)=DSORT(APP) 

APP=RR+SS/VAR3’^=^2 + SSS/VAR3=J^*4 
AP3( J)^=DS0RT( APP) 

APP=RR+SS/VAR4^-v2 + SSS/VAR4*4-4 
AP4( J) =OSORT( APP) 

APP=RR + SS/VAR5~*2+SSS/VAR5=p^4 
AP5( J)=DSORT( APP> 

X=X+.l 

297 CONTINUE 

WRITE(6,1100) (X3(J) ,F( l,Jlt API ( J) ,AP2(J ) , AP3 ( J ) , AP4( J ) 
1,AP5( J) ,0=1,31) 

1100 FORMAT! • »,7F12.7) 

DC 18 K=l,31 
F5(K ) = F( 1,K) 

F6(K)=F( 2,K) 

F7(K)=F (3,K) 

F8(K)=F(4,K) 

F9(K)=FI5,K) 

18 FIO(K) =F (6,K) 

CALL DRAW! 3 1 , X3 , F 5 » 1 , 0 , LABE L , I Tl TLE ,0.5,1.C,0,0,2,2,6, 

16,1, LAST) 

CALL DRAW! 31 , X3 , F6 , 2 , 0 , LABEL , IT I T LE , 0 .5 ,1.0, 0,0, 2, 2, 6, 

16, 1, LAST ) 

CALL DRAW!31 ,X3, F7,2, 0, LABEL, IT I TLE, 0.5, 1.0, 0,0, 2, 2, 6, 

16, 1,LAST) 

CALL DRAW!31,X3, F8 , 2 , 0 „LAB EL , I T I TLE , 0. 5 ,1.0, 0,0, 2, 2, 6, 

16. 1, LAST) 

CALL DRAW! 31,X3,F9,2, 0, LABEL , I TI TLE ,0. 5 , 1 . 0 ,0 ,0 ,2 ,2 ,6, 

16.1, LAST ) 

CALL DRAW! 31 ,X3 , FI 0 , 3 , 0 , LABL5 , I T IT L E , 0 . 5 , 1 . 0 , 0 , 0 , 2 , 2 , 6 

1.6.1, LAST) - 
55 CONTINUE 

STOP 
. END 



C DEGRADATION FACTOR. RC FILTERING AND SAMPLING 

IMPLICIT REALM'S! A-H,Q-Z) 

REAL*8 ITITLEC12)/' S WETT • , 1 l=f= ' */ 

REAL=^4 X3,F5,F6,F7, F3, F9,F10 
REAL*4 LABEL/4H / 

REAL*8 TITLE! 10) /10*» •/ 

DIMENSION Z!43) ,F2 !43 ) ,W I!43 ),Y!43), DELY!43) ,B! 21 ), 
1F!9,31).F5!31),F6!31),F7!31),F8!31),F9!31) ,F10I31 ) , 
2AP1!31),AP2!31 ),AP3!31),AP4!31) ,AP5! 311 ,AP6!31) , 
3F12131) ,F13!31 ) , SB ( 2 1 ) , X3 ! 3 1 ) 

4,FC!9,31),E!20),0!20) 

GAUSS! E, VAR) =0 .3989423^ !OEXP!-! E/ VAR) **2/2. ) ) /VAR 
DO 98 1=1,22 

98 READ! 5,99)WI ! I ) 

99 FCRMAT!F7.1) 

DO 55 MM=1,6 



READ(5,88)N 

88 FCPMAT(I2) 

WRITE( 6,89) N 

89 FORHAK* *,I2) 

RE AD I 5, 999 ) VARl , VA R2 , V AR3 , VAR4 , VARS 
999 FORMAK 5F4. 1) 

CO 5 M=2,6 
DO 1 1=1.22 
READ( 5,2)F2( I ) 

2 FORMAT(F8.5) 

1 CONTINUE 

READ(5,47)A,VAR 

47 FORMAT(2F15.7) 

STP = ( 2/( N-1.0) ) 

FF=0 .0 
DER=0.0 

C CALCULATICN OF DERIVATIVE OF EXPECTED VALUE 

DO 79 1=1. N 

0(1 ) = ( I-(N+1.0 )/2.0)-*=STP 
E( I) = ( I-N/2. 0)*STP + FF 
IF(I.NE.l) GO TO 2000 
E3=E( I ) 

E1=0.0 

E2=GAUSS(E3.VAR) 

GO TO 2500 

2000 IF(I.NE.N) GO TO 2001 
E1=F2 

E2=0.0 
GO TO 2500 

2001 E4=E(I) 

E1=E2 

E2=GAUSS(E4,VAR) 

2500 PPEL=0(I ):^(E1-E2) 

79 DER=D£R+PPEL 

WRITE(6,48 ) A.VAR.DER 

48 FORMAT! • • .3F15.7) 

Z(21)=.98 

Z(22) = l .0 
S=0.0 

CO 3 1=1.20 
Z( I )=S 
S=S+.05 

3 CONTINUE 

C FIT POLYNOMIAL 

CALL ■ LSQPL2(22 .-16.Z.F2.WI . Y . DE LY , B, SB . T IT LE ) 

TRC=20.0 

XX=.01 

DO 35 1=1.31 
NN=TRC/XX-1.0 
RUN=0.0 
RRUN=0.0 
. DO 4 J=1,NN 
T=XX=?=J 
X=DEXP(-T) 

IF(X.LT.O.l) GO TO 3332 

PUN=B( 1)+B(2)*X+E(3 ) *X^*2+B ( 4 ) *X=f^*3 + B ( 5 
1+B( 6)*X^*5 + B(7)^X^=^6+B ( 8)«X*=^ 7+B (9)=J-X^*8+B (10 )=^=X**9 
2+3(11) *X'^” 10+B (12)- X”-^ 1 1+B (13)^ X===’'- 12+3(14)^ X^* 1 3 
3+B(15)^^ X*«14+3(16)*X^«15 + B(17)'-X«=^16 
GO TO 3323 

3332 PUN=F2(3 )*10.0+X 

3333 RUN=RUN+PUN 
RRUN=RRUN+J=5=PUN 

4 CONTINUE 
Al=A*XX/2. 0 
A2=XX*A^RUN 
A3=XX*=?‘2*A*RRUN/TRC 
F41=A1+A2-A3 
F42=A1+A2 

F( M. I ) =DSORT (F42 )/CER 
FC(M.I )=DS0RT(F41) /OER 
IF( I.NE.l) GO TO 66 






XX=XX+.09 
GO TO 35 
66 XX=XX+.l 
35 CONTINUE 
5 r.ONT INUE 
XX = .l 

F(ltl)=1.0 

YYY=1.0-( l.O-DEXP(-TRC) )/TRC 
FC(1,1)=DSQRT(YYY) 

£ 11=09 

C NO QUANTIZER CASE 

DO 7 I=2»31 

YY=( XX/2.0) /DTANH ( XX/2. 0) 

Y4=(XX/2.0)/DSINH(XX/2.0) 

Y3=YY-( 1. D-CEXP (-TRC) )vY4**2/TRC 
F( It I )=DSQRT( YY) 

FC(1 . I )=nSGRT( Y3 ) 

X3( I ) =XX 
XX=XX+.l 
7 CONTINUE 
WRITE! 6,987) 

987 FORMAT ('1') 

WRtTE(6,31) (X3 ( I ) t F(l, I) ,F(2t I) t F(3, I), F(4, I ) ,F(5t I ) , 
1F( 6, I ) , I =1,31) 

X=0.0 

C APPROXIMATION 

DO 297 J=l,31 
RR = F( 1, J )=- *2 

SS=X* (2.0/ (N-1 .0) )=f*2/24.0 

APP=RR+SS/VAR1**2 

APK J)=DSORT( APP ) 

APP=RR+SS/VAR2^-'?2 
AP2( J)=DSORT( APP) 

APP=RR+SS/VAR3^*2 
AP3( J) =DSQRT( APP) 

APP=RR+SS/VAR4-*2 
AP4( J)=OSORT(APP) 

APP=RR+SS/VAR5^*2 
AP5( J)=DSORT( APP) 

X=X+.l 

297 CONTINUE 

WRITE(6, 888) 

888 FORMAT! • •,/////) 

WRITE(6,1100H X3( J) ,F( 1,J) ,AP1( J) ,AP2(J) ,AP3(J) ,AP4(J) 
1 ,AP5( J) , J = l,31 ) 

1100 FORMAT! • • ,7F12.7) 

WRITE!6, 13) 

13 FORMAT! • • ,///// ) 

31 FORMAT! • • ,7F12. 7) 

DO 18 K=l,3i 
F5!K)=F ! 1,K) 

F6!K)=F!2,K) 

F7!K)=F!3,K) 

F8!K)=F!4,K) 

F9!K)=F!5,K) 

18 F10!K)=F!6,K) 

CALL DRAW! 31 ,X3 , F6 , 1 , 0 , LAB EL , I T IT'L E , 0 . 5 , 0 . 0 , 0,0, 1,2, 6, 
16, ItLAST ) 

CALL DRAW!31,X3, F5,2,C,LABEL,ITITLE,0.5,C.C,0,0,1,2,6, 
16, ItLAST) 

CALL DRAW!31,X3, F7 , 2 , 0 , L ABEL , I T I TLE , 0 . 5 , 0. 5 , 0 , 0 , 2 ,2,6, 
16, ItLAST) 

CALL DRAW! 31 , X3 , F 8 , 2 , 0 , LABE L , I T I TLE ,0. 5,0. 5, 0,0, 2, 2, 6, 
16, 1, LAST ) 

CALL DR AW! 31, X3, F9 , 2 , 0 , LABE L , I T IT L E , 0 . 5 , 0 . 5 rO , 0, 2 , 2 , 6, 
16, ItLAST ) 

CALL DRAW!31,X3,F10,3,0,LABEL,ITITLE,0.5,C.5,0,0,2,2,6 
1,6, ItLAST) 

55 CONTINUE 
STOP 
END 
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TRANSFER CHARACTERISTICS OF QUANTIZER WITH OFFSET 
N=NUM6ER OF STEPS 

VAR=NORMALIZED STANDARD DEVIATION OF X 
RHO=CORRELATION COEFICIENT OF X 
DEL=INTEGRATION STEP 
XF=X(N ) 

X0=X(0) 

F=OFFSET 

IMPLICIT REAL«8( A-H»0-Z) 

DIMENSION 0(34) ,X (34) » FUN( 10001 ),W( 18,2) 1 M( 18,2) 

C PROBLEM PARAMETERS 

N=13 

WRITE(6,666)N 
666 FORMAT( • • ,12,//) 

VAR=0.2 

C ALGORITHM PARAMETERS 

RHQ=-0.9999 
D£L=.001 
F=-0.5 
R00=0. 0 
XF=10.0D 10 
X0=-10.3D 10 
SY=-5.0 
SS=0.0 
S=0.0 
B=0.0 

STP=(2./ (N-1.) ) /VAR 
PP=. 3989423 
L=10.0/DEL+1 

C QUANTIZER CHARACTERISTICS 

DO 81 1=1, N 

0( I ) = ( I- (N + 1 .) /2 .0 )*STP 
X(I)=( (I-N/2.0)=S'-STP)+F/VAR 
81 CONTINUE 

C LIMITS OF INTEGRATION 

Nl = l 
N2=N 
N3=N-1 

DO 911 IRR=1,N3 
IF(X( IRR).LE,-5.0) GO TO 1000 
IF(X( IRR) .GE.+5.0 ) GO TO 1001 
GO TO 911 

1000 N1=IRR+1 
GO TO 911 

911 CONTINUE 
GO TO 912 

1001 N2=IRR 

912 DO 91 IR=N1,N2 

IF( IR.NE.Nl ) GO TO 922 
M( IR,1)=1 

K3=(X( IR)+5.0) /DEL+l.O 
.. WO=X(IR)-(K3-1)^DEL+5.0 
M( IR,2)=K3 
W(IR,2)=W0 
GO TO 91 

922 IF( IR.NE.N2) GO TO 92 
M(iR,l)=M( IR-1,2)+1 
W( IR ,1) =CEL-WO 
M( IR,2)=L 
GO TO 91 

92 M(IR,1)=M(IR-1,2)+1 
W( IR ,1 ) = DEL-WO 
K4=(X( IR)+5.0)/DEL+l 
WC=X( IR )-( K4-1 )«DEL+5.0 
M( IR,2) = K4 
W( IR ,2 )=W0 
91 CONTINUE 
C FINAL POINT 

DO 98 1=1, N 
IF( I .NE. 1) GO TO 99 
E2=X0/ 1.414 
E1=X(I ) /1.414 
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GO TO 100 

99 IF( T.NE.N) GO TO 101 
E2=E1 

E1=XF/1.414 
GO TO 100 

101 E2=E1 

E1=X( I ) /1.414 
100 R2 = 0ERF( El )-DERF(E2 ) 

RR=0( I ):<:*2*R2 
SS=SS+RR 
98 CONTINUE 
RATT=SS/2. 0 
RAT=DSORT('RATT ) 

WRITE(6.102)RATT»RAT 

102 FORMAT! 2F15. 9) 

SY=-5.0 

WRITE(6»84)N1,N2 
84 FORMAT! • • » 214) 

DO 30 MM=1,21 
C FORM FUNCTION PHI 

A=2.0'^! 1.0-RH0^!'*2) 

AA=DSORT !A) 

DO 2 J=ltL 
DO 1 I=1»N 
IFII.NE.l) GO TO 5 
E2=! XO-RHO*SY) /AA 
E1 = !X! I )-RHC=PSY)/AA 
GO TO 25 

5 IFII.NE.N) GO TO 3 
E2 = E1 

El = ! XF-RHO^^SY) /AA 
GO TO 25 
3 E2=E1 

El = ! X! I )-RHC*SY) /AA 
25 P,1 = DERF! ED-DERFIE2) 

R=0!I)’^R1 

S=S+R 

1 CONTINUE 

FUN! J) = ! S/ 2.0) *=DEXP! !SY*=^2)/ 1-2.0) ) 
S=.0 

SY=SY+DEL 

2 CONTINUE 

C MAIN INTEGRATION 

DO 7 I1=N1,N2 
IFIIl.NE.Nl) GO TO 111 
M1=M! 11,1) 

M2 = M! 11,2) 

CORR = W! Ii,2)=?'FUN(M2) 

GO TO 15 

111 IF! I1.NE.N2) GO TO 11 
M1=MI 11,1) 

M2=MII1,2) 

CORR=W! II ,1 )*FUN!M1 ) 

GO TO 15 
11 M1 = M!H,1) 

M2=M! 11,2) 

CORR = W! I1,1)=*'FUN!M1) + W! 1 1 ,2 ) *F UN ! M2 ) 
15 M3=M2-1 

DO 21 K0=M1 ,M3 
011=FUNI KO)+FUN! KO+1) 

E=E+011 
21 CONTINUE 

R0=0! I 1 )*! 5+CORR) 

FCQ=R00+B0 
B = .0 

7 CONTINUE 
ROOO=ROO*PP 
POGO-I 1. 0/RATT )>^FOCG 
WRITE! 6,40) RHO,ROOO,ROOO 
40 FORMAT !3C25. 9) 

WRITE! 7,204) RC CO, RCCQ 
204 FORMAT! 2F8. 5) 
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R00=0.0 

S=.0 

SY=-5.0 

IF(MM.NE.l ) GO TO 67 
RHG=RHO+0.0499 
GO TO 30 
67 RhC=RHQ+0.05 
30 CONTINUE 
STOP 
END 



C DEGRAOATICN FACTOR. RC FILTERING AND SAMPLING 

C OFFSET CASE 

IMPLICIT REALT8 ( A-H.O-Z) 

REAL^4 F5.F6.F7.X3 
REAL LA8EL/4H / 

REAL*8 ITITLEI 12 )/ • SWETT •/ 

REAL*8 TITLE! 1 J)/10^' •/ 

DIMENSION Z( 22) ,F2(22) .WK22) »DELVt22) f B(22) ,SB(22) 
l.FC(2.31).APl(31),E(20)t0(20)tF5(31),F6131),F7(31), 
2X3(31) .F(2.31) 

GAUSS (E, VAR )=0.3989A 23* ( DE XP {-( E /VAR) **2/2. ) ) /VAR 
VAR=0.2 ' 

FF=-0. 5 
N=13 

AA=7. 2906748 
G=6. 24225363 
A=AA-G 

STP=(2/(N-1.0) ) 

DER=0.0 

C CALCULATION OF DERIVATIVE OF EXPECTED VALUE 

DO 79 1=1. N 

0(I) = ( I-(N + 1.0)/2.0)*STP 
E( I) = ( I-N/2.0)*STP + FF 
IF(I.NE.l) GO TO 2000 
E3=E(I) 

E1=0 .0 

E2=GAUSS(E3.VAR) 

GO TO 25U0 

2000 IF(I.NE.N) GO TO 2001 
E1=E2 

E2=0.0 
GO TO 2500 

2001 E4=E(I) 

E1=E2 

E2=GAUSS (E4.VAR) 

2500 FPEL=0( I )*( E1-E2) 

79 CER=DER+PPEL 

. WRITE(6.2501)A,DER 

2501 FORMAT!' '.2F12.7) 

CO 98 1=1,21 

98 READ! 5.99) Will) 

99 FORMAT! F 7.1) 

DO 1 1=1.21 
READ! 5,2)F2!I ) 

1 CONTINUE 

2 F0RMATIF8.5) 

S=0.0 

CG 3 1=1.21 
Z!I)=S 
S=S+. 05 

3 CONTINUE 

C FIT POLYNOMIAL 

CALL LS0PL2(21,-11.Z,F2,WI .Y.DELY, B, SB, TITLE) 

TRC=20.0 

B! 1)=0.0 

XX=.01 

CO 35 1=1.31 
NN=TRC/XX-1. 0 
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RUN =0,0 
FRUN=0.0 
DO 4 J = 1,NN 
T=XX=^J 
X=DEXP(-T) 

IF( X.LT. 0. 1) GO TO 3332 

PUN=8( 1 )+B(2)*X + B(3 )=^X*=«=2+S( 4)*X**3+B( 5)*X^^4 
1+B (6)*X^'-i'5 + E (7 )*X^*6+3 (8 J=«X’?:’<=7 + B(9 X^T-S + B ( 10 J'f=X*-T9 
2+B(ll)=i^X=^i'10+B ( l2)TX:i'.icll 
GO TO 3333 

3332 PUN=(F2(3)-F2( 1) )^10.0*X 

3333 RUN=RUN+PUN 
RRUN=RRUN+J*PUN 

4 CONTINUE 
Al=Ai=XX/2.0 
A2=XX^AA*RUN 
A3=XX*=*^2’^AA«RRUN/TRC 
F41=A1+A2-A3 
F42=A1+A2 

FC(2,I )=DSORT( F41) /DER 
F(2. I)=DS0RT(F42)/DER 
IF(I.NE.l) GO TO 66 
XX=XX+. 09 
GO TO 35 

66 XX=XX+.l 
35 CONTINUE 

XX = .l 

F(l,l)=1.0 

YYY=1.0-( l.O-DEXP(-TRC) ) /TRC 
FC(1 ,1 )=CSORT(YYY) 

X3( 1) =0. 0 
CO 7 1=2,31 

YY=(XX/2.0)/DTANH(XX/2 .0) 

Y4=(XX/2.0) /DSINH( XX/2,0) 

Y3=YY-(1 ,0-DEXP (-TRC) ) =?Y4’!'=fc2/TRC 
F(1,I ) =DSORT(YY) 

FC( 1, I )=DSORT( Y3) 

X3( I )=XX 
XX=XX+. 1 
7 CONTINUE 
WRITE(6,987) 

987 FORMAT! • 1* ) 

WRITE(6,31)(X3II),F(l,I),F(2,I),I=lt31) 

X = 0.0 

DO 297 J=l,31 

RR = F(1, J )**2 

SS = (X*STP*-i'2)/24,0 

APP=RR+SS/VAR*^2 

APK J)=DSORT(APP) 

X=X+. 1 

297 CONTINUE 

WRITE(6,888) 

888 FORMAT! • •,/////! 

WRITE!6,31 ) !X3! J ) ,F! 1,J) ,AP1! J) , J = l,31) 

WRITE!6,13) 

13 FORMAT! • • , /////) 

WRITE!6,31)!X3!I),FC!1,I),FC!2,I),I=1,31) 

31 FORMAT! • • ,3F12.7) 

DO 67 K=1,3C 
F5!K ) = F!1,K) 

F6!K)=F! 2,K) 

67 F7!K) = AP1!K) 

CALL DRAW! 30, X 3, F6, 1,0, LAB EL, ITITLE,0.5,O.C,0,0,1,2,6, 

16,1, LAST) 

CALL DRAW! 30,X3, F5, 2, 0,LABEL, ITITLE ,0.5 ,0. C,0,0,1 ,2,6, 

16.1, LAST) 

CALL DRAW! 30, X3,F7, 3, 4, LAB EL, I TITLE, 0.5, 0.0, 0,0,1 ,2,6, 

16.1, LAST ) 

STOP 

END 
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C DYNAMIC RANGE OF I.F. SAMPLING 

IMPLICIT REALMS! A-H,0-Z) 

REAL*4 X3»F5,F6»F7,F8,F9,F10 
REAL LABEL/4H / 

REALMS ITITLEI 12) /' SWETT ',11*' •/ 

DIMENSION 0(34), X(34) ,Y(6,250) ,X3{ 250) ,F5( 250) t 
1F8(250) ,F9(250) , FI 0 (25 0 ) , P ER (25 0 ) , F6 ( 250 ) , F7( 250 ) 
X0=-10.0E 10 
XF=10.0E 10 
00 50 JJ=1 »6 
READ( 5* 99)N 
99 FORMAT (12) 

WRITE(6»88)N 
88 FORMAT( ' • « 12) 

S=.0 
VVAR=. 1 
00 8 L=l,250 
VAR=DSORT( WAR ) 

DO 1 I=1»N 

X( I ) = ( I-N/2. )/VAR 

0(I)=(I-(N+l.)/2.)/VAR 

IF(I.NE.l) GO TO 5 

E2=X0/ 1 .414 

E1=X( I ) /I. 414 

GO TO 6 

5 IF( I .NE.N) GO TO 3 
E2=E1 

E1=XF/1.414 
GO TO 6 
3 E2=E1 

E1=X( I ) /I. 414 

6 R1=DERF( El )-DERF(E2 ) 

R=Q(I )^=«=2^^R1 

S=S+R 

1 CONTINUE 
0ARR=VVARVS 
Y( JJtL )=OARR 
X3 (L)=VVAR 

PER(L) =( 1.0-OARR/(2.0>J-VVAR) )=C^100. 

VVAR=VVAR+.l 

S=.0 

8 CONTINUE 

WRITE(6. 82) ( X3(L) , Y(JJ,L) ,PER(L) »_L = 1 »2 50,3) 

82 FORMATI ' * .3F15.9) 

50 CONTINUE 

DO 46 K=l,250 
F5(K)=Y(1.K) 

F6(K)=Y( 2.K) 

F7(K)=Y(3.K) 

F8(K)=Y(4,K) 

F9(K) = Y( 5,K) 

46 F10(K)=Y(6.K) 

CALL DRAW! 25 0, X3 , F 5 . 1 , 0 ♦ LABEL , I T IT LE »3 . 0 , 5 . 0 » 0» 0 t 2 1 2 , 6 

1.7. 1, LAST) 

CALL DRAW(250,X3,F6,2,0,LABEL, ITITLE,3.0,5.0,0,0,2,2,6 

1.7.1, LA ST) 

CALL DRAW! 250, X3,F7, 2,0, LA BEL, I TITLE, 3. 0,5. 0,0,0, 2, 2, 6 

1.7.1, LAST) 

CALL DRAW! 250, X3 ,F 8 ,2 , 0 , LABE L , I TI TLE, 3. 0,5. 0,0, 0,2, 2, 6 

1,7,1, LAST) 

CALL DRA'rtI 250,X3,F9,2,0,LABEL,ITITLE,3.0,5.0, 0,0,2,2,6 

1.7.1, LAST) 

CALL DR AW ( 2 50, X3, F 10, 3, 0, L AB EL, I TITLE, 3. C, 5. 0,0, 0,2, 2, 

16. 7.1, LAST) 

STOP 

END 
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C DYNAMIC RANGE OF RC FILTERING AND THEN SAMPLING 

IMPLICIT REALMS! A-H,0-Z) 

REAL=^=8 ITITLEI12 )/• SWETT ’/ 

REAL LABEL/4H / 

REAL*4 X3, F5,F6,F7,F8,F9,F10 

DIMENSION 0(34) ,X(34) »Y (7,250) ,X3 (250 ) , F5 ( 250 ) , 
1F8(250),F9(250),F10(250),PER(250) ,NN(8) ,VAN(5) 
2,F6(250) ,F7(250) 

S0 = 5.0 
XF=10.0E 10 
XC=-10.0E 10 

READ(5,73M VAN( I ) ,1=1,1) 

73 FORMAT (F4.1) 

READ(5,71) (NN( I) ,1=1,7) 

71 F0RMAT(7I2) 

1 = 1 

WRITE(6,76)VAN( I ) 

76 FORMAT( ♦ ' ,F4. 1,//) 

VAR=VAN( I) 

DO 70 JJ=1,6 
N=NN( JJ ) 

WRITE(6,72)N 

72 FORMAT( * ♦ ,12,//) 

ALPHA=1.0 

DC 19 K=l,110 
S=0.0 

DO 20 11=1, N 
X(I I) = ( II-N/2.0) 

0( II )^( I I-(N+1 .0)/2.0) 

IFdl.NE.l ) GO TC 99 

E2 = (X0+( VAR-ALPHA=i'VAR)*SO)/ (1.414*VAR*ALPHA) 

E1=(X( I I )+(VAR-ALPHA«VAR)*SQ)/( 1.414* VAR* ALPHA) 

GO TO 100 

99 IF(II.NE.N) GO TO 101 
E2 = E1 

E1=(XF+ (VAR-ALPHA*VAR)*SQ)/ ( 1.414*VARvALPHA J 
GO TO 100 
101 E2=E1 

E1=(X(I I )+ (VAR-ALPHA*VAR)*SO)/ ( 1 .414*VAR* ALPHA ) 

100 R1=DERF(E1)-DERF(E2) 

R2 = 0( I I )*R1 
S=S+R2 
20 CONTINUE 

Y(JJ,K)=S/2.0 
X3(K) =ALPHA 
IF(K.EO.l) GO TO 33 

PER(K)=( 1.0-(Y (JJ,K)/VAR)/( ( ALPHA-1.Q)*S0) )*IC0.0 
GO TO 34 

33 PER(1)=0.0 • 

34 ALPHA=ALPHA+.l 
19 CONTINUE 

WRITE(6,75)(X3(L),Y(JJ,L),PER(L) ,L=1,110,3) 

75 FORMAT (♦ » ,3F15.7) 

70 CONTINUE 
9 CONTINUE 

DO 46 K=l,110 
F5(K) = Y( 1,K) 

F6(K)=Y(2,K) 

F7(K)=Y(3,K) 

F8(K ) = Y( 4,K) 

F9(K)=Y(5,K) 

46 FIOI K)=Y(6,K) 

CALL DRANf 110, X3,F 5, 1,0, LABEL, I TITLE ,2. 0,2. 0,0, 0,2 ,2,6 
1 ,6,1, LAST) 

CALL DRAW! 1 1 0 , X3 , F 6 ,2 , 0 , LA BEl , I T I TL £ , 2 . 0 , 2 .5 ,0,0, 2, 2, 6 

1,7,1, LAST) 

CALL DRAW(110,X3 ,F7,2 ,0,LABEL, IT IT L E , 3 . 0 , 5 . 0, 0, 0 , 2, 2, 6 

1,7,1, LAST) 

CALL DRAW(110,X3,F8,2,0,LABEL,ITITLE,3.0,5.0,0,0,2,2,6 

1,7,1, LAST) 

CALL DRAWI 110, X3,F 9, 2,0, LAB EL, I TITLE, 3. 0,5. 0,0, 0,2, 2, 6 

1,7,1 ,LAST) 
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CALL DRAW( 110,X3 ,F 1 3 ,3 , 0 , LAB EL , IT IT LE ,3 . 0 ♦ 5 . 0 t 0 , 0 , 2 , 2 , 
16,7, 1, LAST) 

STOP 

END 



C DYNAMIC RANG OF N.A.R. RADIOMETER 

IMPLICIT REAL*8( A-H,0-Z) 

REAL=^=4 X3,F5,F6,F7 ,F8,F9,F10 
REAL LABEL/^H / 

REAL*8 ITITLE(12)/' SWETT ',11*' '/ 

DIMENSION 0(34) ,X(34) ,X3(250i ,F5(250) ,F6(250) ,F7(250) 
X0=-10.0E IQ 
XF=10.0E 13 
READ(5,99)N 
99 FORMAT! 12) 

WRITE(6,88)N 
88 FORMAT! • *,I2) 

S = .0 

VVAR=. 1 

DO 8 L=l,250 

VAR=DSORT(VVAR) 

DO 1 1 = 1, N 

X! I ) = ! I-N/2. )/VAR 

0(1) =( I-!N+1.) /2.)/VAR 

IF(I.NE.l) GO TO 5 

E2=X0/ 1.414 

E1=X(I ) /1.414 

GO TO 6 

5 IF(I.NE.N) GO TO 3 
E2=E1 

E1=XF/1.414 
GO TO 6 
3 E2=E1 

E1=X(I ) /I. 414 

6 R1 = DERF! El )-DERF(E2 ) 

R=0(i )^^2*R1 

C +D 

1 CONTINUE 
OARR=VVAR*S 
F5(L ) = 0ARR 
X3!L)=VVAR 
VVAR=VVAR+. 1 
S=.0 

8 CONTINUE 

DO 87 K=l,180 
F6(K ) = F5!K + 40) 

87 F7!K)=F6(K)/F5 (K) 

WRITE!6,82)!X3(L),F5(L),F6!L),F7!L),L=1,180,3) 

82, FORMAT!' ',4 FI 5. 9) 

CALL DRAW!133,X3,F5,1,3,LAB£L,ITITLE,3.3,3.0,3,0,2,2,6 

1.7., 1, LAST) 

CALL DRAW! 180, X3 , F6 ,2 , 0 , LABEL , ITITLE,3.0,5.C,0,0,2,2,6 

1.7.1, LAST) 

CALL DRAW! 180 , X3 , F7 ,3 , 0, LABEL , I TI TL E , 3. 0 , 5 . 0 , 0 , 0 , 2 , 2 ,6 
1 ,7,1 .LAST) 

STOP 

END 
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